Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1161 1316 (N21W14 L=194)   579    225        1 (1051/1290-1257/1344) -30.3 W07-W20 N18-N25
      2 1379  907 (S07W19 L=199)   250     68        1 (1362/ 871-1394/ 944) -90.0 W17-W20 S09-S04
      3 1657  880 (S15W39 L=218)   428    124        0 (1621/ 833-1687/ 931)  38.7 W35-W41 S16-S12
      4 1833  890 (S19W56 L=236)   728     71        2 (1804/ 859-1872/ 924)  31.0 W52-W61 S20-S16
      5 1783  828 (S22W50 L=230)   362     75        1 (1752/ 793-1801/ 860) -68.2 W46-W52 S24-S19
      6 1772  650 (S34W51 L=231)   704    111        1 (1736/ 627-1810/ 682)  -5.7 W46-W56 S36-S31
      7 1518  576 (S32W26 L=205)  1408    267        0 (1455/ 479-1594/ 681)  47.5 W20-W32 S37-S25
      8 1023  623 (S17E08 L=172)   568    191        0 ( 956/ 593-1096/ 659)  -6.3 E11-E02 S19-S14
      9  518  514 (S17E45 L=135)  1460    342        1 ( 370/ 483- 666/ 543)  -4.6 E57-E32 S19-S13
     10  973  362 (S35E18 L=162)  3091    501        7 ( 832/ 204-1130/ 531)  29.6 E39-E02 S46-S24
  

    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