Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  744 1709 (N59E10 L=359)  1605    197        4 ( 676/1645- 822/1751)   0.0 E18-W00 N53-N62
      2 1209 1704 (N46W34 L= 42)  1080    170        1 (1128/1675-1295/1727)  -7.1 W25-W43 N43-N48
      3 1442 1679 (N38W52 L= 60)   766    109        1 (1395/1654-1501/1695) -10.3 W46-W58 N35-N39
      4 1655 1423 (N17W53 L= 62)   459    119        0 (1604/1379-1706/1458) -49.4 W49-W56 N13-N20
      5  605 1339 (N33E24 L=345)   294     86        1 ( 571/1314- 646/1355)   0.0 E26-E20 N31-N33
      6 1334 1257 (N15W24 L= 32)  1037    233        3 (1275/1154-1367/1345)  80.5 W21-W25 N08-N20
      7  758 1219 (N23E14 L=355)   308     98        1 ( 747/1157- 772/1268)  81.9 E14-E12 N18-N25
      8 1879 1094 (S08W65 L= 74)   597     74        1 (1862/1060-1895/1134)  45.0 W62-W68 S09-S05
      9  589 1080 (N17E27 L=342)   482    150        0 ( 521/1051- 650/1112)  12.5 E31-E22 N15-N17
     10 1746  718 (S27W47 L= 56)  2457    245        4 (1687/ 609-1780/ 818) -66.4 W41-W53 S35-S19
     11  933  219 (S48E30 L=339)   527     41        1 ( 905/ 201- 960/ 241)   0.0 E32-E26 S49-S45
  

    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