Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1134 1875 (N57E31 L=212)  1552    145        1 (1074/1832-1190/1918) -10.3 E43-E21 N53-N60
      2  746 1703 (N29E41 L=202)  1747    411        2 ( 657/1581- 845/1834) -20.9 E62-E28 N22-N36
      3 1249 1798 (N54E13 L=230)   380     33        1 (1223/1785-1275/1813)  53.1 E15-E10 N51-N56
      4 1337 1791 (N57E05 L=238)   425     56        0 (1314/1767-1367/1806)  66.8 E07-E02 N53-N60
      5 1160 1557 (N31E07 L=236)   790    137        2 (1110/1513-1210/1608) -23.2 E10-E03 N26-N33
      6 1112 1314 (N14E03 L=240)   400    124        0 (1085/1261-1153/1368)  90.0 E03-E01 N10-N18
      7 1699 1165 (N23W40 L=283)   359    102        0 (1669/1125-1733/1208) -23.2 W36-W43 N21-N24
      8  413 1051 (S18E39 L=204)   527    101        1 ( 357/1031- 455/1066)  14.0 E43-E35 S19-S16
      9  504  968 (S21E31 L=212)   350     54        0 ( 490/ 927- 521/1002) -78.7 E32-E29 S23-S18
     10  830  845 (S19E06 L=237)   582    201        2 ( 765/ 770- 896/ 872)  53.1 E09-E02 S25-S17
     11 1311  725 (S13W25 L=268)   397    127        1 (1270/ 667-1336/ 774) -39.8 W21-W28 S15-S10
  

    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