Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1333 1753 (N57W19 L=200)   467     93        1 (1299/1721-1359/1776)  20.6 W15-W23 N53-N59
      2  935 1318 (N18E09 L=171)   747    305        1 ( 900/1182- 967/1455) -74.5 E12-E05 N09-N27
      3 1792 1388 (N33W57 L=237)   702    194        1 (1733/1356-1827/1410)  -5.7 W50-W60 N30-N34
      4 1865 1131 (N16W63 L=243)   510    136        2 (1842/1088-1885/1193) -74.1 W59-W65 N13-N20
      5 1753  928 (N03W52 L=232)   390     73        1 (1723/ 901-1776/ 954) -31.0 W48-W54 N01-N04
      6 1806  626 (S15W67 L=247)  6165    577       12 (1713/ 514-1884/ 785)  60.9 W58-W76 S22-S04
      7  610  723 (S23E24 L=156)   257     87        1 ( 580/ 682- 652/ 762) -38.7 E26-E20 S24-S20
      8 1755  708 (S10W57 L=238)   577    136        3 (1706/ 674-1783/ 751)  33.7 W53-W60 S12-S07
      9  198  696 (S30E62 L=118)  1040    268        2 ( 181/ 648- 226/ 733)  45.0 E65-E58 S34-S27
     10 1293  636 (S20W22 L=203)   234     95        1 (1247/ 620-1345/ 650)  18.4 W19-W25 S21-S18
     11 1567  536 (S23W46 L=226)   905    173        1 (1497/ 472-1644/ 607)  55.0 W41-W49 S28-S17
     12 1432  274 (S45W50 L=230)   742    112        2 (1382/ 258-1475/ 299)  20.6 W45-W53 S46-S42
  

    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