Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  453 1734 (N25E70 L=217)   888    100        0 ( 426/1697- 500/1767) -71.6 E74-E66 N21-N27
      2  977 1551 (N25E17 L=270)  2079    499        4 ( 873/1318-1047/1717) -49.0 E29-E09 N11-N34
      3  586 1519 (N13E41 L=247)   298     39        1 ( 567/1495- 606/1545) -33.7 E42-E38 N12-N15
      4 1807 1189 (N27W50 L=338)  3922    586        2 (1704/1053-1939/1344) -37.4 W38-W72 N19-N34
      5 1527  960 (N05W31 L=318)   832    220        0 (1461/ 893-1565/1069)  80.5 W28-W32 N00-N12
      6  342 1020 (S20E44 L=244)   488     96        0 ( 321/ 978- 385/1063) -31.0 E46-E39 S22-S17
      7  125  979 (S27E71 L=216)   630     59        1 (  92/ 966- 155/ 997)   7.6 E80-E65 S27-S25
      8 1900  903 (N14W67 L=355)  1106    129        1 (1879/ 853-1927/ 956) -63.4 W62-W72 N10-N16
      9  965  626 (S28W07 L=295)   478    119        0 ( 908/ 596-1026/ 650)  45.0 W04-W10 S31-S25
     10  410  553 (S49E36 L=251)  2071    271        3 ( 310/ 480- 533/ 633)  -8.7 E48-E22 S53-S44
     11 1652  528 (S14W57 L=345)   932    137        0 (1588/ 507-1725/ 555)  31.0 W52-W63 S16-S10
     12 1335  202 (S43W63 L=351)  1655    142        1 (1281/ 161-1390/ 259)   6.0 W53-W78 S45-S39
  

    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