Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1560 1637 (N25W60 L=228)  1211    189        2 (1512/1593-1610/1679)   3.8 W53-W68 N22-N27
      2 1721 1551 (N15W70 L=238)   980    138        1 (1662/1504-1787/1585) -38.7 W64-W76 N10-N18
      3 1305 1474 (N25W31 L=199)   533    178        2 (1218/1447-1387/1491) -31.0 W25-W35 N21-N27
      4 1478 1372 (N14W37 L=206)   291    111        0 (1428/1354-1531/1399) -45.0 W35-W40 N11-N16
      5  312 1219 (N35E48 L=120)   536    164        1 ( 275/1162- 367/1251)  11.3 E52-E41 N31-N36
      6  309 1035 (N23E49 L=119)   592    147        1 ( 269/ 986- 364/1096) -45.0 E52-E44 N19-N27
      7  251  754 (N07E61 L=107)  1285    187        5 ( 201/ 668- 285/ 815)  45.0 E67-E56 N02-N10
      8  462  653 (S02E45 L=123)   367     78        0 ( 437/ 629- 495/ 684) -90.0 E46-E42 S03-N00
      9 1774  587 (S44W49 L=217)   534     53        1 (1748/ 568-1795/ 616) -36.9 W45-W54 S46-S40
     10  380  483 (S11E61 L=107)   628    123        1 ( 344/ 447- 428/ 512)   9.5 E67-E55 S13-S09
  

    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