Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1104 1820 (N46E25 L=328)   386    119        0 (1066/1785-1134/1851) -16.7 E30-E20 N43-N47
      2 1477 1623 (N44W12 L=  5)   753    123        1 (1433/1583-1515/1651)  29.7 W08-W15 N41-N47
      3 1547 1606 (N46W19 L= 12)   387     62        1 (1526/1576-1569/1633)  90.0 W16-W20 N43-N49
      4 1710 1383 (N36W35 L= 29)  3688    475        4 (1587/1209-1785/1536) -40.1 W23-W43 N25-N44
      5  396 1484 (N04E54 L=299)   477    115        0 ( 360/1445- 420/1508)  -9.5 E57-E51 N01-N05
      6 1049 1262 (N07E05 L=348)   228     84        0 (1035/1221-1069/1303) -53.1 E06-E02 N04-N09
      7  172 1261 (S13E68 L=285)  1050    202        2 ( 128/1224- 214/1295)  10.3 E74-E63 S14-S10
      8 1552 1162 (N16W27 L= 20)   436    121        0 (1506/1102-1578/1226) -35.5 W22-W29 N13-N18
      9 1806 1100 (N22W48 L= 41)   448    102        0 (1776/1079-1837/1133)  -9.5 W44-W51 N20-N23
     10 1470 1012 (N05W25 L= 19)   819    206        2 (1405/ 964-1562/1087)  -4.4 W19-W32 N03-N07
     11 1873  970 (N17W58 L= 52)   788    104        2 (1852/ 921-1889/1025) -63.4 W55-W61 N14-N20
     12 1714  938 (N09W44 L= 38)  1189    262        2 (1620/ 910-1801/ 960)  21.8 W37-W52 N04-N11
     13 1725  756 (N00W52 L= 45)   929    165        1 (1657/ 696-1777/ 831)  73.3 W48-W54 S05-N05
     14  737  749 (S30E09 L=344)   197     53        0 ( 711/ 729- 759/ 776)  63.4 E10-E08 S31-S27
     15 1068  686 (S24W13 L=  6)   547    106        1 (1001/ 659-1117/ 714)   7.1 W07-W15 S25-S21
     16 1505  578 (S16W43 L= 36)   711    200        0 (1412/ 512-1589/ 632)  63.4 W39-W46 S22-S10
     17 1214  578 (S26W25 L= 19)   485    122        1 (1169/ 532-1259/ 624) -12.5 W20-W29 S27-S23
  

    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