Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  854 1487 (N30E05 L=159)   261     88        0 ( 819/1446- 893/1511) -51.3 E07-E02 N27-N32
      2 1427 1471 (N22W33 L=197)   367     68        1 (1402/1431-1452/1507) -90.0 W31-W34 N18-N24
      3  243 1214 (N21E54 L=109)  1107    167        2 ( 190/1156- 316/1257)  12.1 E60-E46 N17-N23
      4 1726  975 (S12W46 L=210)   462    157        2 (1676/ 951-1788/1014) -15.3 W41-W52 S13-S09
      5 1405  836 (S16W21 L=185)  1337    406        1 (1353/ 663-1446/ 996)  78.7 W18-W23 S26-S06
      6 1804  899 (S17W54 L=217)   938    122        3 (1767/ 867-1855/ 927) -11.3 W49-W59 S19-S15
      7 1747  852 (S20W48 L=212)   469     40        1 (1717/ 828-1774/ 875) -31.0 W45-W50 S21-S17
      8 1811  806 (S23W55 L=219)   361     69        1 (1789/ 770-1830/ 841) -63.4 W52-W58 S25-S21
      9 1297  752 (S20W13 L=177)   270    120        0 (1259/ 705-1316/ 805)  50.2 W10-W15 S22-S16
     10  266  545 (S19E68 L= 95)  1519    158        5 ( 238/ 482- 297/ 627)  39.3 E76-E61 S23-S14
     11 1676  419 (S50W55 L=219)   826     80        0 (1636/ 403-1720/ 440)  -9.5 W47-W65 S51-S47
  

    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