Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  961 1791 (N60W19 L= 13)   404     53        1 ( 933/1767- 987/1819) -45.0 W16-W23 N57-N62
      2  853 1744 (N59W04 L=358)  1553    222        4 ( 784/1668- 927/1794)  25.2 E03-W15 N53-N62
      3 1322 1756 (N46W49 L= 43)  1484    184        3 (1230/1722-1404/1785) -11.9 W39-W59 N43-N48
      4 1447 1742 (N41W60 L= 54)   447     54        1 (1425/1714-1474/1765)  63.4 W56-W62 N39-N43
      5 1519 1731 (N38W66 L= 61)   831     85        2 (1476/1707-1556/1750) -10.3 W60-W73 N36-N39
      6 1633 1571 (N26W61 L= 55)   411     60        1 (1606/1536-1654/1598)  76.0 W59-W63 N23-N27
      7 1519 1335 (N15W38 L= 33)  1392    281        2 (1432/1237-1565/1431) -81.9 W35-W41 N08-N22
      8  785 1381 (N33E09 L=345)   258     74        1 ( 755/1359- 824/1396)   0.0 E11-E06 N31-N33
      9  773 1129 (N17E14 L=340)   736    303        0 ( 636/1093- 896/1173)  -3.2 E23-E05 N15-N18
     10  253 1090 (N21E56 L=298)   565     71        0 ( 226/1064- 288/1124) -29.7 E59-E51 N19-N23
     11 1602  699 (S25W33 L= 27)   369     86        0 (1577/ 677-1641/ 730) -38.7 W31-W36 S26-S22
     12 1125  276 (S46E10 L=344)   990    174        1 (1048/ 242-1194/ 315)  12.1 E17-E03 S48-S43
  

    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