Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1716 1460 (N27W57 L= 22)   393    132        1 (1671/1447-1755/1472) -14.0 W52-W60 N25-N27
      2 1136 1386 (N27W09 L=335)  1035    336        2 ( 958/1332-1272/1429) -22.8 E02-W18 N22-N30
      3 1404 1269 (N17W26 L=352)  1065    241        6 (1334/1206-1477/1332) -45.0 W22-W30 N13-N21
      4 1751 1261 (N14W54 L= 20)   532    110        0 (1723/1209-1782/1310) -56.3 W52-W56 N10-N16
      5  680 1181 (N15E22 L=304)   408    178        2 ( 644/1091- 717/1264) -70.0 E24-E19 N09-N20
      6 1781  885 (S09W54 L= 19)  1714    345        2 (1694/ 788-1857/ 954) -35.2 W45-W63 S16-S04
      7  451  903 (S01E38 L=287)   272     47        2 ( 427/ 886- 472/ 921)  33.7 E40-E36 S02-N00
      8 1423  816 (S11W25 L=350)   362     77        1 (1401/ 770-1447/ 863)  71.6 W23-W26 S13-S07
      9 1516  575 (S27W33 L=359)  2323    498        4 (1332/ 427-1666/ 754)  45.0 W20-W44 S37-S16
     10 1227  623 (S22W11 L=337)   424    128        1 (1163/ 575-1295/ 671)  29.1 W07-W16 S24-S19
  

    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