Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1524 1687 (N27W61 L=259)   741     77        3 (1489/1658-1557/1718)  68.2 W56-W64 N24-N29
      2 1286 1617 (N31W38 L=236)  1691    379        1 (1146/1569-1457/1672) -43.2 W29-W47 N24-N39
      3  975 1593 (N40W16 L=214)   447    102        0 ( 924/1569-1020/1614)  -7.1 W11-W19 N38-N41
      4 1459 1445 (N17W39 L=238)   821    296        2 (1296/1425-1562/1479) -21.3 W28-W46 N12-N20
      5 1636 1373 (N07W48 L=247)   515     71        0 (1609/1344-1662/1409) -78.7 W46-W50 N04-N09
      6 1054 1330 (N22W10 L=209)   650    181        1 ( 970/1305-1134/1367) -15.3 W05-W16 N20-N24
      7 1529 1296 (N07W38 L=236)   275     68        0 (1494/1276-1563/1317) -45.0 W35-W39 N04-N08
      8 1108 1240 (N16W11 L=209)   311    165        0 (1041/1213-1182/1258) -14.0 W06-W14 N13-N16
      9  505  993 (N17E32 L=166)   737    286        1 ( 399/ 959- 645/1017) -22.4 E40-E23 N11-N19
     10  344  940 (N17E46 L=152)   279     22        0 ( 325/ 915- 361/ 959)  63.4 E47-E44 N16-N18
     11 1748  906 (S22W43 L=242)   339     75        0 (1708/ 891-1786/ 920) -20.6 W39-W47 S23-S20
     12  327  777 (N08E52 L=146)   961    164        3 ( 297/ 724- 376/ 842) -66.8 E55-E48 N05-N12
     13 1460  786 (S20W19 L=217)   292     97        0 (1435/ 743-1485/ 825)  38.7 W16-W21 S21-S17
     14 1357  696 (S22W10 L=209)   248     95        0 (1339/ 653-1376/ 745)  45.0 W08-W12 S23-S18
     15 1104  640 (S18E06 L=192)   383    127        3 (1065/ 594-1138/ 703)  74.1 E08-E04 S20-S13
     16 1553  417 (S50W21 L=219)   419     84        0 (1516/ 391-1579/ 450)   8.1 W16-W23 S51-S47
     17 1356  304 (S51E01 L=197)  1074    184        0 (1291/ 268-1437/ 345)  -8.1 E07-W07 S53-S48
     18  481  307 (S24E68 L=130)   621     36        1 ( 461/ 285- 502/ 331)  38.7 E72-E63 S25-S21
  

    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