Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1152 1819 (N48W45 L= 31)   376     53        0 (1134/1789-1164/1846)  35.5 W42-W49 N45-N50
      2 1544 1695 (N27W62 L= 48)   826     96        1 (1486/1674-1590/1710) -29.1 W56-W67 N24-N29
      3 1342 1578 (N27W38 L= 23)   555    146        2 (1283/1548-1413/1613) -54.5 W35-W41 N23-N30
      4 1625 1507 (N14W54 L= 40)  1481    307        1 (1499/1454-1751/1540) -31.4 W44-W64 N07-N19
      5 1135 1403 (N23W17 L=  3)   842    227        2 (1025/1358-1238/1437) -39.3 W11-W22 N17-N27
      6 1623 1293 (N03W44 L= 29)   296     73        1 (1590/1260-1651/1320)  21.8 W41-W46 N01-N04
      7  171 1006 (N23E64 L=282)  1310    152        2 ( 107/ 975- 254/1032)  -2.5 E75-E52 N21-N25
      8 1384  999 (S06W19 L=  5)   215     50        0 (1356/ 978-1404/1025)  33.7 W17-W21 S06-S04
      9 1300  842 (S12W11 L=356)   232     49        1 (1274/ 818-1324/ 867)   0.0 W09-W12 S13-S11
     10 1484  627 (S31W18 L=  3)  1749    292        6 (1433/ 508-1535/ 736)  52.4 W12-W23 S38-S23
     11 1607  543 (S41W27 L= 13)   590     68        2 (1578/ 515-1643/ 574)   0.0 W24-W31 S43-S38
     12  925  323 (S34E29 L=316)   417     77        0 ( 910/ 286- 943/ 365)  50.2 E31-E26 S37-S31
  

    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