Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  757 1723 (N58W07 L=343)  2961    497        4 ( 595/1586- 946/1860)   8.4 E14-W40 N52-N62
      2 1575 1673 (N24W65 L= 41)  6530    414        8 (1431/1594-1711/1758) -50.5 W55-W79 N16-N34
      3 1669 1532 (N14W61 L= 38)  1310    243        0 (1547/1474-1752/1576) -42.7 W53-W68 N08-N20
      4 1323 1502 (N24W34 L= 10)   240     79        0 (1285/1490-1359/1525)  -9.5 W30-W37 N22-N25
      5  719 1309 (N30E11 L=325)   310     82        1 ( 674/1285- 757/1338)  -9.5 E14-E08 N29-N31
      6  950 1212 (N18W01 L=337)   218     71        1 ( 931/1183- 970/1251)  45.0 W00-W03 N16-N20
      7 1828  828 (S30W52 L= 28)  1189     95        2 (1777/ 803-1876/ 858) -23.2 W46-W60 S33-S26
      8 1402  735 (S21W14 L=351)   469    143        0 (1385/ 664-1414/ 805)  63.4 W12-W16 S25-S16
      9 1207  585 (S24E02 L=335)   558    220        1 (1154/ 504-1254/ 680)  42.0 E07-W03 S28-S19
     10 1756  602 (S44W45 L= 21)   521     44        0 (1731/ 576-1780/ 623)   9.5 W41-W48 S45-S42
     11 1341  527 (S33W05 L=341)   300     60        1 (1324/ 490-1360/ 554)  36.9 W03-W07 S34-S30
     12 1339  394 (S43W01 L=337)   722    141        2 (1301/ 360-1397/ 452) -78.7 E01-W05 S47-S37
     13  906  398 (S28E27 L=309)   539    114        1 ( 863/ 364- 964/ 429) -66.8 E28-E24 S31-S24
     14  533  316 (S24E61 L=275)  1415     93        1 ( 496/ 274- 576/ 355) -78.7 E66-E56 S27-S20
     15  689  290 (S30E49 L=287)   600     72        2 ( 669/ 259- 708/ 330)  45.0 E53-E45 S32-S26
     16  701  236 (S34E54 L=282)   785    102        2 ( 652/ 216- 738/ 258)  -5.7 E59-E48 S36-S32
  

    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