; You should have all the binary individual images representing each region from the
; "vornoi-bse2-pass1.pro" file.
; For this file I have all those images names as "Region_XX.tif"

;PROCEDURE regiongrow START
pro regiongrow ,x,y,im

IF im[x,y] eq 254 THEN BEGIN
      im[x,y] = 255
      regiongrow ,x-1,y-1,im
      regiongrow ,x,y-1,im
      regiongrow ,x+1,y-1,im
      regiongrow ,x-1,y,im
      regiongrow ,x+1,y,im
      regiongrow ,x-1,y+1,im
      regiongrow ,x,y+1,im
      regiongrow ,x+1,y+1,im
END

END
;PROCEDURE regiongrow END

;PROCEDURE distancemap START
pro distancemap,image,dist

dist = MORPH_DISTANCE(image,/BACKGROUND, NEIGHBOR_SAMPLING = 3);
    for i=0,767 do begin
       for j=0,1023 do begin
         IF dist[j,i] EQ 0 THEN BEGIN
          dist[j,i]=-1
         END
       end
    end

end
;PROCEDURE distancemap END

pro bse2_distance_transform

device , decomposed = 0
loadct , 0

file1 = filepath('Region_1.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image1 = read_image(file1)           ;THIS IS TO READ A IMAGE FILE.

file2 = filepath('Region_2.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image2 = read_image(file2)           ;THIS IS TO READ A IMAGE FILE.

file3 = filepath('Region_3.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image3 = read_image(file3)           ;THIS IS TO READ A IMAGE FILE.

file4 = filepath('Region_4.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image4 = read_image(file4)           ;THIS IS TO READ A IMAGE FILE.

file5 = filepath('Region_5.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image5 = read_image(file5)           ;THIS IS TO READ A IMAGE FILE.

file6 = filepath('Region_6.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image6 = read_image(file6)           ;THIS IS TO READ A IMAGE FILE.

file7 = filepath('Region_7.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image7 = read_image(file7)           ;THIS IS TO READ A IMAGE FILE.

file8 = filepath('Region_8.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image8 = read_image(file8)           ;THIS IS TO READ A IMAGE FILE.

file9 = filepath('Region_9.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image9 = read_image(file9)           ;THIS IS TO READ A IMAGE FILE.

file10 = filepath('Region_10.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image10 = read_image(file10)           ;THIS IS TO READ A IMAGE FILE.

file11 = filepath('Region_11.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image11 = read_image(file11)           ;THIS IS TO READ A IMAGE FILE.

file12 = filepath('Region_12.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image12 = read_image(file12)           ;THIS IS TO READ A IMAGE FILE.

file13 = filepath('Region_13.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image13 = read_image(file13)           ;THIS IS TO READ A IMAGE FILE.

file14 = filepath('Region_14.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image14 = read_image(file14)           ;THIS IS TO READ A IMAGE FILE.

file15 = filepath('Region_15.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image15 = read_image(file15)           ;THIS IS TO READ A IMAGE FILE.

file16 = filepath('Region_16.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image16 = read_image(file16)           ;THIS IS TO READ A IMAGE FILE.

file17 = filepath('Region_17.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image17 = read_image(file17)           ;THIS IS TO READ A IMAGE FILE.

file18 = filepath('Region_18.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image18 = read_image(file18)           ;THIS IS TO READ A IMAGE FILE.

file19 = filepath('Region_19.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image19 = read_image(file19)           ;THIS IS TO READ A IMAGE FILE.

file20 = filepath('Region_20.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image20 = read_image(file20)           ;THIS IS TO READ A IMAGE FILE.

file21 = filepath('Region_21.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image21 = read_image(file21)           ;THIS IS TO READ A IMAGE FILE.

file22 = filepath('Region_22.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image22 = read_image(file22)           ;THIS IS TO READ A IMAGE FILE.

file23 = filepath('Region_23.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image23 = read_image(file23)           ;THIS IS TO READ A IMAGE FILE.

file24 = filepath('Region_24.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image24 = read_image(file24)           ;THIS IS TO READ A IMAGE FILE.

file25 = filepath('Region_25.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image25 = read_image(file25)           ;THIS IS TO READ A IMAGE FILE.

file26 = filepath('Region_26.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image26 = read_image(file26)           ;THIS IS TO READ A IMAGE FILE.

file27 = filepath('Region_27.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image27 = read_image(file27)           ;THIS IS TO READ A IMAGE FILE.

file28 = filepath('Region_28.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image28 = read_image(file28)           ;THIS IS TO READ A IMAGE FILE.

file29 = filepath('Region_29.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image29 = read_image(file29)           ;THIS IS TO READ A IMAGE FILE.

file30 = filepath('Region_30.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image30 = read_image(file30)           ;THIS IS TO READ A IMAGE FILE.

file31 = filepath('Region_31.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image31 = read_image(file31)           ;THIS IS TO READ A IMAGE FILE.

file32 = filepath('Region_32.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image32 = read_image(file32)           ;THIS IS TO READ A IMAGE FILE.

file33 = filepath('Region_33.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image33 = read_image(file33)           ;THIS IS TO READ A IMAGE FILE.

file34 = filepath('Region_34.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image34 = read_image(file34)           ;THIS IS TO READ A IMAGE FILE.

file35 = filepath('Region_35.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image35 = read_image(file35)           ;THIS IS TO READ A IMAGE FILE.

file36 = filepath('Region_36.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image36 = read_image(file36)           ;THIS IS TO READ A IMAGE FILE.

file37 = filepath('Region_37.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image37 = read_image(file37)           ;THIS IS TO READ A IMAGE FILE.

file38 = filepath('Region_38.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image38 = read_image(file38)           ;THIS IS TO READ A IMAGE FILE.

file39 = filepath('Region_39.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image39 = read_image(file39)           ;THIS IS TO READ A IMAGE FILE.

file40 = filepath('Region_40.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image40 = read_image(file40)           ;THIS IS TO READ A IMAGE FILE.

file41 = filepath('Region_41.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image41 = read_image(file41)           ;THIS IS TO READ A IMAGE FILE.

file42 = filepath('Region_42.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image42 = read_image(file42)           ;THIS IS TO READ A IMAGE FILE.

file43 = filepath('Region_43.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image43 = read_image(file43)           ;THIS IS TO READ A IMAGE FILE.

file44 = filepath('Region_44.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image44 = read_image(file44)           ;THIS IS TO READ A IMAGE FILE.

file45 = filepath('Region_45.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image45 = read_image(file45)           ;THIS IS TO READ A IMAGE FILE.

file46 = filepath('Region_46.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image46 = read_image(file46)           ;THIS IS TO READ A IMAGE FILE.

file47 = filepath('Region_47.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image47 = read_image(file47)           ;THIS IS TO READ A IMAGE FILE.

file48 = filepath('Region_48.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image48 = read_image(file48)           ;THIS IS TO READ A IMAGE FILE.

file49 = filepath('Region_49.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image49 = read_image(file49)           ;THIS IS TO READ A IMAGE FILE.

file50 = filepath('Region_50.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image50 = read_image(file50)           ;THIS IS TO READ A IMAGE FILE.

file51 = filepath('Region_51.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image51 = read_image(file51)           ;THIS IS TO READ A IMAGE FILE.

file52 = filepath('Region_52.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image52 = read_image(file52)           ;THIS IS TO READ A IMAGE FILE.

file53 = filepath('Region_53.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image53 = read_image(file53)           ;THIS IS TO READ A IMAGE FILE.

file54 = filepath('Region_54.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image54 = read_image(file54)           ;THIS IS TO READ A IMAGE FILE.

file55 = filepath('Region_55.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image55 = read_image(file55)           ;THIS IS TO READ A IMAGE FILE.

file56 = filepath('Region_56.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image56 = read_image(file56)           ;THIS IS TO READ A IMAGE FILE.

file57 = filepath('Region_57.tif' , $
subdirectory = ['idldata1','bse2','bse2-vernoi','Regions'])
imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image57 = read_image(file57)           ;THIS IS TO READ A IMAGE FILE.

cnt=0
dist=[1024,768]
distancemap,image1,dist
dist1=dist

distancemap,image2,dist
dist2=dist

distancemap,image3,dist
dist3=dist
;
distancemap,image4,dist
dist4=dist

distancemap,image5,dist
dist5=dist

distancemap,image6,dist
dist6=dist

distancemap,image7,dist
dist7=dist

distancemap,image8,dist
dist8=dist

distancemap,image9,dist
dist9=dist

distancemap,image10,dist
dist10=dist

distancemap,image11,dist
dist11=dist

distancemap,image12,dist
dist12=dist

distancemap,image13,dist
dist13=dist

distancemap,image14,dist
dist14=dist

distancemap,image15,dist
dist15=dist

distancemap,image16,dist
dist16=dist

distancemap,image17,dist
dist17=dist

distancemap,image18,dist
dist18=dist

distancemap,image19,dist
dist19=dist

distancemap,image20,dist
dist20=dist
;
distancemap,image21,dist
dist21=dist

distancemap,image22,dist
dist22=dist

distancemap,image23,dist
dist23=dist

distancemap,image24,dist
dist24=dist

distancemap,image25,dist
dist25=dist

distancemap,image26,dist
dist26=dist

distancemap,image27,dist
dist27=dist

distancemap,image28,dist
dist28=dist

distancemap,image29,dist
dist29=dist

distancemap,image30,dist
dist30=dist

distancemap,image31,dist
dist31=dist

distancemap,image32,dist
dist32=dist

distancemap,image33,dist
dist33=dist

distancemap,image34,dist
dist34=dist

distancemap,image35,dist
dist35=dist

distancemap,image36,dist
dist36=dist

distancemap,image37,dist
dist37=dist

distancemap,image38,dist
dist38=dist

distancemap,image39,dist
dist39=dist

distancemap,image40,dist
dist40=dist

distancemap,image41,dist
dist41=dist

distancemap,image42,dist
dist42=dist

distancemap,image43,dist
dist43=dist

distancemap,image44,dist
dist44=dist

distancemap,image45,dist
dist45=dist

distancemap,image46,dist
dist46=dist

distancemap,image47,dist
dist47=dist

distancemap,image48,dist
dist48=dist

distancemap,image49,dist
dist49=dist

distancemap,image50,dist
dist50=dist

distancemap,image51,dist
dist51=dist

distancemap,image52,dist
dist52=dist

distancemap,image53,dist
dist53=dist

distancemap,image54,dist
dist54=dist

distancemap,image55,dist
dist55=dist

distancemap,image56,dist
dist56=dist

distancemap,image57,dist
dist57=dist


;dist1 = MORPH_DISTANCE(image1,/BACKGROUND, NEIGHBOR_SAMPLING = 3);
;
;for i=0,767 do begin
;    for j=0,1023 do begin
;       IF dist1[j,i] EQ 0 THEN BEGIN
;            dist1[j,i]=-1
;       END
;    end
;end

minarray = intarr(57,1)
index = indgen(1024,768)

for i=0,767 do begin
    for j=0,1023 do begin
       minarray[0]=dist1[j,i]
       minarray[1]=dist2[j,i]
       minarray[2]=dist3[j,i]
       minarray[3]=dist4[j,i]
       minarray[4]=dist5[j,i]
       minarray[5]=dist6[j,i]
       minarray[6]=dist7[j,i]
       minarray[7]=dist8[j,i]
       minarray[8]=dist9[j,i]
       minarray[9]=dist10[j,i]
       minarray[10]=dist11[j,i]
       minarray[11]=dist12[j,i]
       minarray[12]=dist13[j,i]
       minarray[13]=dist14[j,i]
       minarray[14]=dist15[j,i]
       minarray[15]=dist16[j,i]
       minarray[16]=dist17[j,i]
       minarray[17]=dist18[j,i]
       minarray[18]=dist19[j,i]
       minarray[19]=dist20[j,i]
       minarray[20]=dist21[j,i]
       minarray[21]=dist22[j,i]
       minarray[22]=dist23[j,i]
       minarray[23]=dist24[j,i]
       minarray[24]=dist25[j,i]
       minarray[25]=dist26[j,i]
       minarray[26]=dist27[j,i]
       minarray[27]=dist28[j,i]
       minarray[28]=dist29[j,i]
       minarray[29]=dist30[j,i]
       minarray[30]=dist31[j,i]
       minarray[31]=dist32[j,i]
       minarray[32]=dist33[j,i]
       minarray[33]=dist34[j,i]
       minarray[34]=dist35[j,i]
       minarray[35]=dist36[j,i]
       minarray[36]=dist37[j,i]
       minarray[37]=dist38[j,i]
       minarray[38]=dist39[j,i]
       minarray[39]=dist40[j,i]
       minarray[40]=dist41[j,i]
       minarray[41]=dist42[j,i]
       minarray[42]=dist43[j,i]
       minarray[43]=dist44[j,i]
       minarray[44]=dist45[j,i]
       minarray[45]=dist46[j,i]
       minarray[46]=dist47[j,i]
       minarray[47]=dist48[j,i]
       minarray[48]=dist49[j,i]
       minarray[49]=dist50[j,i]
       minarray[50]=dist51[j,i]
       minarray[51]=dist52[j,i]
       minarray[52]=dist53[j,i]
       minarray[53]=dist54[j,i]
       minarray[54]=dist55[j,i]
       minarray[55]=dist56[j,i]
       minarray[56]=dist57[j,i]

       result = min(minarray)
       for k=0,56 do begin
         if minarray[k] eq result then begin
            index[j,i]=k+1
         end
       end
    end
end

; Assign different GrayScale value to each vornoi region so that all regions are
; seperated by a distinct boundary.

for i=0,767 do begin
    for j=0,1023 do begin
        if index[j,i] EQ 1 then begin
               image1[j,i] = 5
        end

        if index[j,i] EQ 2 then begin
               image1[j,i] = 55
        end

        if index[j,i] EQ 3 then begin
               image1[j,i] = 105
        end

        if index[j,i] EQ 4 then begin
               image1[j,i] = 205
        end

        if index[j,i] EQ 5 then begin
               image1[j,i] = 10
        end

        if index[j,i] EQ 6 then begin
               image1[j,i] = 60
        end

        if index[j,i] EQ 7 then begin
               image1[j,i] = 110
        end

        if index[j,i] EQ 8 then begin
               image1[j,i] = 210
        end

        if index[j,i] EQ 9 then begin
               image1[j,i] = 15
        end

        if index[j,i] EQ 10 then begin
               image1[j,i] = 65
        end

        if index[j,i] EQ 11 then begin
               image1[j,i] = 115
        end

        if index[j,i] EQ 12 then begin
               image1[j,i] = 215
        end

        if index[j,i] EQ 13 then begin
               image1[j,i] = 20
        end

        if index[j,i] EQ 14 then begin
               image1[j,i] = 70
        end

        if index[j,i] EQ 15 then begin
               image1[j,i] = 120
        end

        if index[j,i] EQ 16 then begin
               image1[j,i] = 220
        end

        if index[j,i] EQ 17 then begin
               image1[j,i] = 25
        end

        if index[j,i] EQ 18 then begin
               image1[j,i] = 75
        end

        if index[j,i] EQ 19 then begin
               image1[j,i] = 125
        end

        if index[j,i] EQ 20 then begin
               image1[j,i] = 225
        end

        if index[j,i] EQ 21 then begin
               image1[j,i] = 30
        end

        if index[j,i] EQ 22 then begin
               image1[j,i] = 80
        end

        if index[j,i] EQ 23 then begin
               image1[j,i] = 130
        end

        if index[j,i] EQ 24 then begin
               image1[j,i] = 230
        end

        if index[j,i] EQ 25 then begin
               image1[j,i] = 35
        end

        if index[j,i] EQ 26 then begin
               image1[j,i] = 85
        end

        if index[j,i] EQ 27 then begin
               image1[j,i] = 135
        end

        if index[j,i] EQ 28 then begin
               image1[j,i] = 235
        end

        if index[j,i] EQ 29 then begin
               image1[j,i] = 40
        end

        if index[j,i] EQ 30 then begin
               image1[j,i] = 90
        end

        if index[j,i] EQ 31 then begin
               image1[j,i] = 140
        end

        if index[j,i] EQ 32 then begin
               image1[j,i] = 240
        end

        if index[j,i] EQ 33 then begin
               image1[j,i] = 45
        end

        if index[j,i] EQ 34 then begin
               image1[j,i] = 95
        end

        if index[j,i] EQ 35 then begin
               image1[j,i] = 145
        end

        if index[j,i] EQ 36 then begin
               image1[j,i] = 245
        end

       if index[j,i] EQ 37 then begin
               image1[j,i] = 50
        end

       if index[j,i] EQ 38 then begin
               image1[j,i] = 100
        end

       if index[j,i] EQ 39 then begin
               image1[j,i] = 155
        end

       if index[j,i] EQ 40 then begin
               image1[j,i] = 175
        end

       if index[j,i] EQ 41 then begin
               image1[j,i] = 160
        end

       if index[j,i] EQ 42 then begin
               image1[j,i] = 180
        end

       if index[j,i] EQ 43 then begin
               image1[j,i] = 165
        end

       if index[j,i] EQ 44 then begin
               image1[j,i] = 185
        end

       if index[j,i] EQ 45 then begin
               image1[j,i] = 170
        end

       if index[j,i] EQ 46 then begin
               image1[j,i] = 190
        end

       if index[j,i] EQ 47 then begin
               image1[j,i] = 155
        end

       if index[j,i] EQ 48 then begin
               image1[j,i] = 175
        end

       if index[j,i] EQ 49 then begin
               image1[j,i] = 160
        end

       if index[j,i] EQ 50 then begin
               image1[j,i] = 180
        end

       if index[j,i] EQ 51 then begin
               image1[j,i] = 165
        end

       if index[j,i] EQ 52 then begin
               image1[j,i] = 185
        end

       if index[j,i] EQ 53 then begin
               image1[j,i] = 170
        end

       if index[j,i] EQ 54 then begin
               image1[j,i] = 190
        end

       if index[j,i] EQ 55 then begin
               image1[j,i] = 155
        end

       if index[j,i] EQ 56 then begin
               image1[j,i] = 175
        end

       if index[j,i] EQ 57 then begin
               image1[j,i] = 160
        end

          ;   PRINT, j,i,dist50[j,i]
    end
end

WINDOW, 0, XSIZE=imagesize[0], YSIZE=imagesize[1], $
TITLE='ouput'
im1 = rotate(image1,7)
TV , im1
write_tiff, 'vornoi_bse2.tif' , image1        ;TO WRITE FINAL OUTPUT IMAGE.

;file1 = filepath('output.tif' , $
;subdirectory = ['idldata1','bse2','bse2-vernoi'])
;imagesize =[1024,768]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
;origimage = read_image(file1)           ;THIS IS TO READ A IMAGE FILE.
;origimage1=origimage
;
;image3=image1
;image2=image1
;
;for i=0,767 do begin
;    for j=0,1023 do begin
;       image2[j,i]=100
;    end
;end
;
;greencolor = !d.table_size-1
;tvlct ,0,255,0,greencolor
;;image1 = bytscl(image1,top=greencolor-1)
;
;for i=0,767 do begin
;    for j=0,1023 do begin
;       if origimage[j,i] eq 200 then begin
;         image2[j,i]=greencolor
;       end
;    end
;end
;
;redcolor = !d.table_size-2
;tvlct ,255,0,0,redcolor
;
;;image1 = bytscl(image1,top=redcolor-1)
;
;for i=0,767 do begin
;    for j=0,1023 do begin
;       if origimage[j,i] eq 50 then begin
;         image2[j,i]=redcolor
;       end
;    end
;end
;
;bluecolor = !d.table_size-1
;tvlct ,0,0,255,bluecolor
;
;;image1 = bytscl(image1,top=bluecolor-1)
;
;for i=0,767 do begin
;    for j=0,1023 do begin
;       if origimage[j,i] eq 100 then begin
;         image2[j,i]=bluecolor
;       end
;    end
;end
;
;for i=0,767 do begin
;    for j=0,1023 do begin
;       IF ((image2[j,i] EQ redcolor)||(image2[j,i] EQ greencolor)||(image2[j,i] EQ bluecolor)) THEN BEGIN
;           image1[j,i]= image2[j,i]
;        END
;    END
;END

;void = TVRead(/tiff, File='myimage')

END