Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1357 1583 (N41W10 L=299)   265     36        1 (1338/1556-1372/1604) -56.3 W08-W11 N39-N42
      2 1294 1473 (N31W07 L=296)   526    108        1 (1243/1440-1332/1500)   8.1 W03-W10 N29-N32
      3 1196 1329 (N20W04 L=293)   252     82        2 (1151/1312-1232/1345)  31.0 W01-W06 N17-N20
      4 1762 1286 (N31W49 L=338)   549     83        1 (1723/1256-1786/1316)  -8.1 W44-W51 N28-N32
      5 1894 1092 (N22W66 L=355)   745    124        0 (1872/1042-1928/1127)  50.2 W62-W73 N18-N24
      6 1320 1094 (N08W16 L=305)   207     49        0 (1299/1065-1340/1121) -33.7 W14-W17 N06-N09
      7  259 1032 (S17E54 L=235)   330     51        2 ( 229/1010- 283/1056)  -9.5 E57-E51 S17-S15
      9 1745  910 (N07W50 L=339)   356    118        0 (1687/ 890-1778/ 934)  18.4 W44-W53 N05-N09
     10  403  680 (S36E38 L=250)   879    166        3 ( 333/ 643- 472/ 738)  -8.1 E45-E31 S38-S33
     11 1430  498 (S24W41 L=330)   909    264        4 (1326/ 444-1535/ 545)  36.9 W35-W48 S29-S20
  

    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