Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1007 1794 (N52W28 L= 20)   510     60        0 ( 969/1778-1042/1813) -45.0 W24-W31 N49-N54
      2 1567 1683 (N27W62 L= 54)   683     94        0 (1521/1658-1608/1707) -33.7 W59-W66 N24-N29
      3 1500 1599 (N24W49 L= 41)   587    175        0 (1416/1581-1573/1615) -28.6 W43-W54 N21-N27
      4 1586 1506 (N16W50 L= 42)   298    125        0 (1562/1450-1602/1561)  32.0 W45-W53 N14-N19
      5 1701 1489 (N12W60 L= 52)  1381    151        4 (1638/1423-1732/1539) -69.4 W55-W63 N08-N16
      6 1013 1513 (N32W12 L=  4)   704    198        0 ( 921/1487-1114/1535) -26.6 W06-W18 N27-N35
      7 1367 1501 (N22W34 L= 26)   460    147        0 (1287/1487-1436/1514) -16.7 W28-W38 N19-N23
      8  437 1069 (N19E36 L=316)   377     73        2 ( 407/1035- 461/1096)  45.0 E37-E33 N17-N20
      9 1708  642 (S37W40 L= 32)   395     52        2 (1678/ 622-1733/ 670) -15.9 W36-W43 S38-S35
     10 1506  573 (S36W21 L= 13)   371     67        1 (1474/ 547-1538/ 599)  21.8 W18-W23 S37-S34
     11 1371  453 (S40W08 L=  0)   433     77        0 (1330/ 427-1412/ 478)   0.0 W04-W11 S41-S38
     12  599  277 (S32E57 L=295)  1034     97        1 ( 565/ 248- 642/ 316) -90.0 E61-E53 S34-S28
  

    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