Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1602 1636 (N21W60 L= 49)   958    157        0 (1531/1593-1668/1673) -52.1 W56-W63 N16-N25
      2 1428 1613 (N26W45 L= 35)   362     78        0 (1404/1569-1460/1645)  78.7 W43-W47 N23-N28
      3 1302 1570 (N28W35 L= 24)   631    129        0 (1232/1545-1366/1592) -32.0 W30-W38 N25-N30
      4  400 1109 (N27E38 L=312)   464    132        1 ( 363/1063- 457/1141)  12.5 E41-E32 N24-N28
      5  405  994 (N20E39 L=310)   406     74        0 ( 368/ 968- 439/1014) -45.0 E41-E36 N17-N21
      6 1748  879 (S23W42 L= 31)  2307    438        2 (1645/ 817-1860/ 982)   2.7 W33-W54 S26-S16
      7 1154  698 (S15E02 L=348)   457    131        3 (1115/ 644-1189/ 736)  26.6 E05-W01 S17-S12
      8 1334  573 (S28W06 L=355)   838    254        1 (1236/ 485-1421/ 691)  26.6 E02-W14 S31-S23
  

    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