Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  442 1709 (N20E65 L=131)  1830    276        1 ( 418/1660- 490/1757) -45.0 E71-E61 N15-N23
      2  688 1014 (S15E19 L=177)   612    129        0 ( 623/ 977- 748/1060)  54.5 E22-E16 S18-S11
      3 1829  837 (N09W59 L=255)  1419    177        5 (1796/ 773-1864/ 897) -26.6 W54-W66 N06-N12
      4  814  835 (S22E07 L=189)   269     54        0 ( 787/ 809- 846/ 868) -11.3 E09-E04 S23-S20
      5 1412  761 (S09W30 L=226)   679    177        0 (1388/ 674-1430/ 834) -45.0 W25-W33 S13-S05
      6  990  763 (S21W06 L=202)   322    113        2 ( 933/ 742-1055/ 796)   7.1 W01-W09 S22-S20
      7 1340  712 (S14W27 L=223)   375    142        0 (1259/ 675-1385/ 741)   0.0 W21-W30 S15-S12
      8  587  640 (S41E18 L=179)   837    219        2 ( 488/ 622- 685/ 670)  13.2 E26-E09 S43-S37
      9 1423  458 (S25W44 L=240)  2552    922        5 (1102/ 406-1718/ 532)  15.5 W20-W74 S32-S15
     10 1220  256 (S43W46 L=242)   420     55        0 (1190/ 236-1248/ 272)  11.3 W43-W48 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