Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1322 1748 (N40W51 L=241)   484     62        0 (1284/1734-1355/1761) -15.9 W48-W55 N38-N41
      2 1005 1668 (N46W22 L=211)   261     48        0 ( 974/1656-1038/1678) -26.6 W19-W24 N44-N47
      3 1565 1593 (N23W57 L=246)   687    108        0 (1511/1556-1605/1633) -71.6 W54-W58 N19-N26
      4 1508 1469 (N18W44 L=234)  1743    464        1 (1334/1319-1671/1551) -59.4 W35-W50 N05-N27
      5  154  965 (N21E69 L=120)  1031    134        1 ( 121/ 911- 203/1021) -28.3 E75-E62 N17-N24
      6 1907  924 (S25W65 L=255)   852     61        1 (1881/ 894-1936/ 949)   6.3 W60-W72 S26-S23
      7 1326  853 (S10W13 L=203)   276    110        0 (1301/ 806-1350/ 903) -81.9 W12-W13 S13-S06
      8 1642  535 (S42W32 L=222)  1549    221        1 (1550/ 472-1713/ 613)  10.0 W22-W39 S44-S37
      9  469  422 (S17E56 L=133)   788    176        1 ( 439/ 349- 515/ 499)  65.6 E59-E52 S22-S11
     10  736  264 (S34E46 L=143)   490     56        0 ( 703/ 245- 769/ 285) -36.9 E48-E43 S35-S31
  

    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