Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  457 1709 (N20E63 L=131)  1723    276        1 ( 433/1658- 512/1761) -38.7 E70-E59 N16-N24
      2  711 1003 (S16E17 L=177)   520    131        0 ( 644/ 966- 759/1043)  54.5 E20-E14 S19-S12
      3 1840  830 (N09W60 L=255)  1322    184        4 (1807/ 769-1866/ 892) -36.9 W55-W66 N06-N12
      4  833  825 (S22E05 L=189)   192     35        0 ( 812/ 803- 854/ 852) -14.0 E07-E03 S23-S21
      5 1434  769 (S08W31 L=225)   685    101        1 (1414/ 710-1454/ 836) -50.2 W28-W33 S11-S04
      6 1029  745 (S21W08 L=203)   307    131        0 ( 964/ 728-1104/ 773)  11.3 W03-W13 S21-S19
      7  614  629 (S40E15 L=179)   783    196        0 ( 518/ 605- 710/ 652)  14.0 E23-E07 S43-S37
      8 1184  481 (S31W28 L=223)   279    115        0 (1137/ 448-1234/ 513)   0.0 W23-W32 S32-S30
      9 1422  451 (S25W44 L=239)  1052    276        0 (1284/ 427-1529/ 472)  33.7 W36-W52 S30-S20
     10 1229  251 (S43W47 L=241)   398     81        0 (1192/ 230-1262/ 267)  18.4 W43-W50 S44-S42
  

    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