Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1517 1490 (N21W46 L=185)   294    115        0 (1461/1475-1570/1506) -20.6 W41-W49 N18-N22
      2 1117 1425 (N28W16 L=156)   325     90        1 (1091/1391-1141/1463)  36.9 W14-W18 N25-N30
      3 1536 1395 (N15W42 L=182)   251     67        1 (1508/1377-1570/1419) -53.1 W40-W44 N13-N17
      4  142 1026 (N22E71 L= 68)   798    106        0 ( 117/ 977- 159/1073) -63.4 E77-E68 N19-N25
      5 1462  739 (S20W20 L=160)   332    103        0 (1433/ 695-1491/ 782)  31.0 W17-W22 S22-S17
      6 1717  646 (S35W41 L=181)   456     38        1 (1697/ 618-1734/ 669) -63.4 W39-W43 S36-S32
      7  657  411 (S23E40 L= 99)  2370    489        4 ( 500/ 275- 781/ 595)  23.6 E64-E24 S31-S12
      8  456  428 (S18E57 L= 83)  1316    278        3 ( 421/ 314- 489/ 547)  46.7 E66-E49 S27-S10
      9  944  389 (S30E22 L=118)   312     80        0 ( 901/ 371- 976/ 414)   8.1 E25-E18 S31-S29
  

    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