Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1145 1727 (N53W11 L=211)   421     76        1 (1113/1696-1187/1745)   7.1 W07-W15 N49-N54
      2  514 1529 (N35E43 L=157)  1306    266        1 ( 414/1487- 626/1592)  -8.7 E57-E31 N31-N38
      3  720 1474 (N32E23 L=177)   308     97        0 ( 681/1446- 752/1512) -29.7 E27-E20 N29-N34
      4  979 1146 (N11E03 L=197)   194    109        1 ( 960/1098-1005/1191) -71.6 E04-E01 N07-N13
      5  776  685 (S18E16 L=184)   468    213        1 ( 679/ 664- 873/ 720)  -8.7 E22-E09 S19-S16
      6  416  568 (S28E47 L=153)  2201    340        3 ( 331/ 412- 498/ 685)  28.8 E77-E36 S41-S19
      7  830  608 (S23E13 L=187)   348    178        0 ( 779/ 546- 904/ 651)  37.9 E16-E07 S27-S20
      8 1593  423 (S37W52 L=252)   620     83        1 (1571/ 386-1612/ 459) -36.9 W46-W56 S40-S34
      9 1477  362 (S42W42 L=242)  1171    124        2 (1417/ 335-1541/ 390)  14.0 W36-W49 S44-S39
     10 1249  310 (S46W22 L=222)  1170    223        3 (1158/ 288-1336/ 343)   3.4 W13-W30 S48-S43
  

    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