Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  459 1616 (N22E53 L= 27)   726     90        2 ( 438/1582- 494/1665) -36.9 E57-E48 N19-N25
      2  827 1631 (N28E25 L= 55)   243     44        0 ( 815/1599- 839/1663) -59.0 E26-E22 N25-N30
      3 1145 1541 (N27E02 L= 78)   791    247        0 (1056/1461-1243/1596) -24.8 E08-W06 N23-N29
      4  499 1510 (N15E44 L= 36)   605     92        1 ( 462/1482- 546/1539)  -7.1 E47-E39 N13-N16
      5 1703 1267 (N21W41 L=120)   871    259        3 (1661/1162-1758/1376) -53.1 W37-W46 N15-N27
      6  521 1286 (N02E36 L= 44)   263     59        0 ( 491/1262- 549/1311) -26.6 E38-E33 N00-N02
      7  146 1221 (S05E71 L=  9)   597     76        0 ( 112/1206- 191/1241)   8.7 E77-E64 S05-S03
      8  506 1120 (S08E34 L= 46)   804    287        1 ( 443/1047- 582/1233)  61.7 E37-E29 S13-N00
      9 1818 1144 (N16W53 L=132)   614    106        1 (1781/1104-1838/1199)  81.9 W48-W55 N13-N20
     10 1719 1056 (N09W44 L=123)   428    104        0 (1677/1032-1760/1090)  39.8 W40-W47 N06-N11
     11  266  962 (S20E53 L= 27)   405     75        0 ( 249/ 921- 280/1003) -68.2 E54-E51 S22-S17
     12  688  929 (S17E19 L= 60)   208     28        0 ( 659/ 917- 715/ 942)  18.4 E20-E17 S17-S16
     13  723  665 (S33E14 L= 65)  1971    528        6 ( 633/ 496- 842/ 903)  69.2 E20-E08 S45-S16
     14  537  766 (S29E29 L= 51)   669    260        0 ( 461/ 701- 584/ 866)  53.1 E34-E25 S34-S22
     15 1475  844 (S09W31 L=110)   212     75        0 (1430/ 833-1514/ 856)  26.6 W27-W33 S10-S07
     16  323  794 (S30E48 L= 32)   282     31        1 ( 306/ 768- 335/ 820) -90.0 E49-E46 S31-S28
     17 1151  675 (S25W14 L= 94)   243     85        1 (1113/ 650-1184/ 705)  45.0 W12-W16 S26-S22
     18 1627  649 (S16W48 L=127)   301     37        0 (1606/ 629-1649/ 669)  71.6 W45-W49 S17-S14
     19  899  616 (S33E01 L= 78)   479    150        0 ( 822/ 601- 970/ 630)  12.5 E06-W03 S34-S30
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y