Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  791 1740 (N59E08 L= 50)  1269    275        3 ( 658/1704- 916/1774)  -1.9 E24-W06 N56-N60
      2 1092 1414 (N29W11 L= 68)   506    190        0 (1031/1337-1135/1468) -68.2 W07-W12 N23-N33
      3  869 1373 (N30E05 L= 53)   416    139        0 ( 820/1341- 933/1415) -40.6 E08-E01 N26-N33
      4 1571 1314 (N15W41 L= 99)   250     63        1 (1548/1288-1598/1333) -33.7 W39-W42 N13-N15
      5  251 1214 (N26E57 L=  1)   466    163        1 ( 216/1143- 284/1279) -38.7 E63-E52 N21-N29
      6 1592 1245 (N10W40 L= 98)   470    125        0 (1545/1219-1650/1268) -20.6 W37-W45 N08-N12
      7  427 1125 (N19E39 L= 19)   246     59        1 ( 412/1097- 446/1163) -76.0 E40-E37 N17-N21
      8  492 1082 (N16E34 L= 24)   250     79        0 ( 449/1070- 534/1096)   0.0 E37-E30 N15-N16
      9 1910  992 (S12W68 L=126)   666     54        0 (1898/ 956-1924/1026) -78.7 W66-W71 S14-S09
     10  192  937 (N09E63 L=355)   412     76        1 ( 171/ 906- 209/ 972)  76.0 E65-E60 N07-N11
     11 1226  743 (S13W08 L= 66)   254     56        2 (1201/ 722-1260/ 764)   0.0 W06-W10 S14-S11
     12 1732  590 (S34W49 L=107)  2176    386        0 (1587/ 471-1869/ 695)  14.0 W35-W69 S40-S29
     13  406  579 (S14E50 L=  8)   840    191        0 ( 356/ 504- 434/ 666)  42.5 E58-E45 S19-S08
     14  256  603 (S12E66 L=352)   654     67        1 ( 224/ 574- 286/ 619) -26.6 E71-E61 S14-S11
     15 1296  370 (S41W10 L= 68)   393     69        0 (1260/ 352-1336/ 387)   0.0 W06-W13 S42-S39
     16 1109  274 (S47E07 L= 50)  1191    140        2 (1052/ 233-1184/ 305)  13.0 E13-W00 S49-S42
     17  937  203 (S50E26 L= 31)  1445    141        3 ( 896/ 149- 986/ 257)  30.5 E36-E19 S56-S44
  

    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