Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1215 1640 (N39W32 L=  2)  1216    181        8 (1122/1592-1281/1690) -59.0 W27-W36 N34-N44
      2 1689 1541 (N20W64 L= 35)   865    116        3 (1654/1489-1728/1591)  74.1 W61-W68 N15-N22
      3 1716 1400 (N12W57 L= 28)   641     86        1 (1685/1360-1744/1440)  18.4 W52-W62 N10-N13
      4  495 1332 (N36E31 L=299)  1041    235        3 ( 378/1289- 590/1363) -16.7 E43-E23 N31-N39
      5  806 1202 (N22E10 L=321)  1067    485        3 ( 674/1042- 964/1305) -55.0 E17-E02 N09-N30
      6  739 1071 (N15E16 L=314)   405    241        0 ( 681/1002- 790/1115) -63.4 E19-E13 N10-N18
      7 1801  972 (S15W52 L= 22)  1249    270        2 (1745/ 861-1850/1077)  32.7 W45-W59 S20-S09
      8 1769 1033 (S10W50 L= 20)   360     40        0 (1745/1009-1790/1053) -36.9 W47-W51 S11-S08
      9  170  946 (N16E67 L=263)  2300    154        4 ( 136/ 878- 212/1011) -63.4 E72-E61 N11-N19
     10 1703  948 (S13W42 L= 12)   366    106        2 (1678/ 910-1726/ 992)  68.2 W40-W44 S15-S10
     11 1289  833 (S10W12 L=342)   256    143        0 (1234/ 798-1326/ 889)  23.2 W07-W14 S11-S06
     12  905  759 (S06E12 L=318)   423    165        2 ( 847/ 721- 958/ 787)   7.1 E16-E08 S08-S04
     13  673  340 (S29E43 L=288)   496    119        1 ( 618/ 326- 732/ 359) -18.4 E47-E38 S31-S27
  

    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