Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1114 1501 (N28W17 L=334)   358     93        0 (1070/1470-1156/1522) -45.0 W14-W19 N25-N30
      2 1826 1369 (N03W65 L= 22)  1376    160        0 (1766/1315-1896/1405)   2.6 W55-W78 N01-N04
      3 1487 1359 (N11W35 L=353)   365     81        2 (1445/1340-1528/1378) -31.0 W33-W38 N09-N12
      4  614 1158 (N19E22 L=295)   220     75        0 ( 580/1136- 644/1188)  11.3 E24-E19 N17-N19
      5 1705 1019 (S13W42 L=359)   464     85        2 (1675/ 979-1742/1054)  31.0 W39-W45 S14-S11
      6 1719  859 (S22W41 L=358)  1347    213        2 (1639/ 796-1784/ 939)  19.7 W34-W48 S24-S18
      7  477  899 (N07E35 L=282)   285     70        0 ( 439/ 875- 508/ 917) -26.6 E37-E33 N05-N08
      8 1222  864 (S10W08 L=325)   407    233        0 (1166/ 836-1309/ 912)   0.0 W04-W13 S11-S07
      9  309  841 (N06E50 L=267)   739    121        2 ( 260/ 796- 360/ 874) -26.6 E54-E46 N03-N08
     10  233  793 (N05E59 L=258)   614     95        0 ( 206/ 754- 272/ 826)  20.6 E62-E53 N03-N06
     11 1310  648 (S25W09 L=326)  2430    734        5 (1054/ 414-1567/ 820)  17.6 E13-W28 S34-S19
     12  925  688 (S14E13 L=304)   198     87        1 ( 911/ 645- 940/ 737)  68.2 E14-E11 S16-S11
     13  778  446 (S26E30 L=287)   271     92        0 ( 751/ 407- 813/ 487) -90.0 E30-E28 S28-S22
     14  706  372 (S29E39 L=279)   307     51        1 ( 683/ 349- 734/ 395)   0.0 E41-E36 S30-S27
     15  557  279 (S32E59 L=258)  1254    113        6 ( 529/ 221- 596/ 317)  17.5 E70-E51 S36-S29
     16  848  244 (S42E38 L=279)   451     75        1 ( 813/ 208- 884/ 272)  18.4 E42-E33 S44-S39
  

    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