Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1310 1729 (N54W19 L= 83)   342     47        0 (1282/1715-1337/1742)   0.0 W15-W22 N52-N55
      2  390 1606 (N32E64 L=360)   848    114        1 ( 368/1567- 417/1653) -22.6 E72-E58 N29-N34
      3  607 1474 (N25E35 L= 29)  1628    380        2 ( 516/1393- 665/1617) -79.4 E41-E30 N19-N35
      4 1507 1402 (N30W31 L= 96)  2462    546        2 (1289/1231-1630/1527) -33.2 W14-W40 N19-N37
      5 1369 1194 (N15W21 L= 85)   368    158        1 (1339/1147-1410/1268) -54.5 W18-W23 N12-N19
      6 1550 1201 (N17W34 L= 98)   376    123        0 (1526/1140-1564/1251) -74.1 W31-W34 N13-N20
      7  453  780 (S19E37 L= 27)   662    158        0 ( 431/ 708- 491/ 843)  66.0 E38-E33 S23-S14
      8 1465  650 (S18W34 L= 98)   867    216        2 (1384/ 577-1538/ 712)  51.3 W29-W38 S23-S13
      9 1721  646 (S16W56 L=120)   394     50        0 (1705/ 620-1740/ 669)  71.6 W54-W58 S17-S14
     10 1727  567 (S21W62 L=126)  1577    174        4 (1667/ 515-1805/ 609)  21.0 W54-W71 S24-S18
  

    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