Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  602 1657 (N19E47 L=133)  1266    172        4 ( 558/1604- 668/1692)  66.8 E51-E42 N15-N22
      2  250 1467 (N01E71 L=110)   579     85        1 ( 227/1440- 267/1490) -14.0 E74-E66 N00-N02
      3 1743 1178 (N25W41 L=221)   501     87        1 (1704/1143-1776/1216)  -8.1 W37-W44 N23-N25
      4  247  919 (S32E51 L=129)   916    156        2 ( 221/ 843- 270/ 989)  73.3 E52-E48 S36-S26
      5 1840  921 (N14W57 L=237)   457     71        0 (1818/ 893-1859/ 947)  63.4 W54-W58 N11-N15
      6  903  892 (S16E03 L=177)   663    180        1 ( 817/ 838- 991/ 944)  59.0 E06-W00 S21-S11
      7 1617  701 (S06W46 L=226)  1192    200        3 (1575/ 616-1647/ 786) -63.4 W42-W49 S10-N00
      8 1182  656 (S21W20 L=200)   210     94        0 (1138/ 637-1223/ 677)   0.0 W17-W23 S22-S20
      9  707  559 (S42E06 L=174)  1101    292        3 ( 579/ 535- 859/ 589)  18.4 E16-W05 S45-S38
     10 1482  368 (S27W54 L=234)   901    248        0 (1376/ 344-1577/ 397)  36.9 W47-W60 S31-S22
  

    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