Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1547 1639 (N43W47 L=131)  2335    400        8 (1468/1575-1598/1695) -60.9 W39-W54 N38-N48
      2 1115 1511 (N35W07 L= 90)   510    168        0 (1093/1437-1148/1563)  77.5 W05-W09 N29-N38
      3 1713 1323 (N21W51 L=134)   454    100        1 (1675/1280-1744/1357) -56.3 W47-W53 N17-N23
      4  765 1312 (N21E17 L= 66)   235     92        0 ( 745/1269- 783/1351) -59.0 E18-E15 N18-N23
      5 1712 1240 (N16W49 L=133)   716    184        1 (1656/1214-1767/1274)  11.3 W44-W54 N14-N17
      6 1103 1186 (N14W05 L= 88)   319    168        1 (1087/1103-1126/1264)  73.3 W03-W06 N08-N18
      7 1759 1103 (N07W52 L=136)   638    241        1 (1673/1067-1826/1154) -14.9 W44-W59 N04-N10
      8  802 1115 (N09E14 L= 70)   299     90        1 ( 758/1074- 843/1147) -38.7 E16-E11 N06-N11
      9 1345 1121 (N09W20 L=104)   218     73        2 (1316/1092-1379/1147) -36.9 W18-W22 N07-N10
     10  682  969 (S00E21 L= 62)   353    134        1 ( 626/ 935- 753/1010)  24.0 E25-E16 S02-N02
     11  246  782 (S13E58 L= 25)   388     70        0 ( 220/ 758- 273/ 804) -21.8 E61-E55 S14-S11
     12 1590  771 (S13W38 L=122)   330     89        1 (1558/ 756-1628/ 791)   9.5 W35-W41 S13-S11
     13 1723  624 (S23W53 L=136)   810    106        1 (1687/ 582-1754/ 674)  56.3 W49-W55 S26-S20
     14  674  586 (S25E24 L= 59)   456    131        0 ( 613/ 566- 734/ 605)  12.5 E28-E19 S26-S23
     15 1235  369 (S41W17 L=101)  1794    356        6 (1152/ 302-1334/ 476)  54.5 W10-W23 S47-S32
     16  887  363 (S42E11 L= 72)   454     86        1 ( 856/ 343- 916/ 397)  51.3 E14-E09 S43-S38
     17  704  324 (S45E29 L= 55)   379    137        1 ( 671/ 294- 739/ 360)  33.7 E33-E24 S48-S42
  

    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