Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1645 1636 (N22W69 L=246)   506     84        1 (1610/1599-1671/1671)  59.0 W66-W71 N19-N24
      2 1493 1594 (N25W51 L=228)  1003    159        0 (1451/1542-1558/1624)   0.0 W44-W56 N23-N27
      3 1682 1495 (N14W60 L=237)  1457    274        0 (1587/1397-1776/1561) -65.0 W55-W65 N05-N20
      4 1192 1407 (N24W21 L=198)   301    114        0 (1135/1388-1236/1430) -45.0 W18-W23 N22-N27
      5 1612 1411 (N12W49 L=226)   304     84        0 (1562/1400-1649/1422) -23.2 W44-W51 N10-N13
      6  230 1028 (N24E58 L=120)  1309    272        1 ( 192/ 941- 293/1161) -60.3 E62-E51 N18-N32
      7 1828  996 (S18W54 L=231)   444     61        1 (1803/ 966-1860/1028)  18.4 W50-W58 S18-S16
      8 1528  953 (S10W28 L=205)   348     65        0 (1506/ 916-1548/ 990) -90.0 W26-W28 S12-S07
      9  197  796 (N10E67 L=111)  3995    230        7 ( 140/ 676- 270/ 900)  70.3 E76-E58 N02-N16
     10 1188  671 (S17W01 L=178)   212     54        1 (1171/ 641-1207/ 702)  71.6 E01-W01 S18-S15
     11 1812  632 (S42W54 L=231)  1001    120        3 (1775/ 611-1852/ 658)  -7.6 W47-W63 S44-S38
     12 1702  545 (S44W39 L=216)   507     38        0 (1679/ 527-1727/ 561)   9.5 W36-W42 S45-S41
  

    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