Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1511 1514 (N18W45 L= 54)  1063    248        0 (1420/1488-1620/1536) -23.2 W39-W53 N14-N21
      2 1808 1408 (N03W67 L= 75)   934    251        0 (1771/1372-1857/1441)  -3.2 W60-W78 N01-N04
      3 1114 1253 (N16W11 L= 19)   202     58        2 (1077/1237-1146/1274) -53.1 W09-W12 N14-N18
      4 1596 1165 (S03W37 L= 45)   685    166        0 (1566/1089-1650/1248)  83.7 W35-W40 S08-N01
      5 1819  967 (S21W50 L= 58)  1272    334        1 (1769/ 873-1903/1078)  17.5 W43-W62 S25-S17
      6  160  966 (N22E63 L=305)   779    138        2 ( 125/ 930- 190/ 988) -45.0 E66-E58 N20-N23
      7 1745  818 (S27W40 L= 48)   295    104        0 (1712/ 789-1776/ 838) -39.8 W37-W43 S30-S25
      8 1525  739 (S24W22 L= 30)   243     76        0 (1514/ 703-1537/ 780)  78.7 W21-W22 S27-S21
      9 1689  716 (S32W34 L= 42)   303     74        0 (1669/ 678-1710/ 751)  36.9 W32-W36 S33-S30
     10  195  716 (N07E69 L=299)  1028    169        3 ( 154/ 678- 238/ 750) -26.6 E73-E63 N04-N09
     11  493  549 (S08E47 L=321)   382     60        1 ( 469/ 519- 516/ 583) -90.0 E47-E45 S09-S05
  

    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