Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1085 1897 (N57E38 L=206)   515     46        0 (1065/1871-1103/1921) -38.7 E44-E34 N53-N61
      2 1154 1864 (N57E28 L=217)   600     59        0 (1125/1836-1178/1885) -16.7 E32-E22 N54-N59
      3  725 1716 (N29E44 L=201)  1768    404        1 ( 637/1590- 822/1851) -21.0 E68-E29 N22-N37
      4 1233 1806 (N54E15 L=230)   363     41        0 (1208/1794-1259/1824)  36.9 E17-E12 N51-N56
      5 1320 1792 (N57E07 L=238)   379     55        0 (1296/1768-1348/1808)  63.4 E08-E04 N52-N59
      6 1142 1578 (N32E09 L=236)   495     98        1 (1112/1530-1169/1616) -39.8 E12-E06 N28-N33
      7 1084 1323 (N14E05 L=240)   320    112        0 (1063/1270-1109/1375) -81.9 E05-E03 N10-N17
      8 1676 1183 (N23W37 L=282)   282     77        0 (1650/1146-1701/1222) -31.0 W34-W39 N21-N24
      9  407 1054 (S18E40 L=205)   539    129        0 ( 346/1022- 463/1075)   0.0 E45-E34 S19-S16
     10  481  935 (S23E32 L=213)   431     77        0 ( 461/ 886- 506/ 977)  80.5 E32-E30 S26-S20
     11  807  863 (S19E08 L=237)   590    141        2 ( 759/ 815- 871/ 888)  45.0 E11-E04 S22-S16
     12  553  864 (S26E25 L=220)   328     88        0 ( 516/ 842- 599/ 886)  38.7 E27-E22 S27-S23
     13 1288  733 (S13W23 L=268)   409    135        1 (1244/ 677-1314/ 785) -39.8 W19-W26 S16-S10
  

    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