Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  637 1793 (N38E59 L=145)   970    127        0 ( 593/1759- 685/1841)  -8.4 E77-E50 N35-N40
      2  600 1679 (N28E50 L=155)  1021    215        0 ( 516/1641- 675/1721)  41.6 E55-E46 N24-N32
      3 1428 1533 (N39W16 L=221)   935    214        0 (1356/1460-1461/1619) -42.0 W08-W19 N33-N44
      4  321 1541 (N14E68 L=136)   626     83        0 ( 298/1504- 346/1584) -78.7 E70-E66 N11-N16
      5  892 1502 (N23E20 L=185)   226    125        0 ( 856/1458- 928/1558) -32.0 E24-E16 N21-N26
      6 1323 1416 (N28W10 L=215)  1246    472        0 (1186/1225-1402/1550) -38.7 E02-W18 N17-N35
      7  927 1350 (N15E13 L=192)   786    274        4 ( 805/1302-1052/1393)   0.0 E21-E04 N12-N16
      8 1260 1330 (N21W07 L=212)   324    142        1 (1227/1268-1277/1380) -63.4 W04-W09 N17-N23
      9 1674 1221 (N25W39 L=244)   882    176        2 (1648/1150-1706/1278) -60.3 W36-W42 N20-N28
     10  296 1117 (S12E52 L=153)  1413    377        2 ( 255/1000- 335/1257) -60.3 E58-E46 S19-S03
     11 1753 1229 (N27W47 L=252)   328     95        0 (1709/1201-1787/1257)  -6.3 W42-W51 N26-N28
     12  206 1125 (S13E62 L=143)  2847    289        8 ( 145/ 991- 257/1230) -70.3 E71-E54 S20-S06
     13 1881 1024 (N18W65 L=270)  1052    154        1 (1867/ 947-1896/1096) -83.7 W62-W68 N13-N22
     14  365 1001 (S17E43 L=162)   534    161        0 ( 334/ 937- 384/1091)  84.3 E45-E41 S21-S11
     15 1293  879 (S05W19 L=224)   207     72        1 (1254/ 862-1330/ 897)  31.0 W16-W21 S06-S03
     16 1097  699 (S20W12 L=217)   612    217        3 ( 999/ 672-1189/ 730)  28.6 W06-W17 S22-S15
     17  437  647 (S39E35 L=170)  3403    489        8 ( 223/ 573- 624/ 708)   4.7 E66-E17 S44-S33
     18  779  388 (S49E00 L=204)   654    113        1 ( 736/ 345- 831/ 434) -24.0 E05-W05 S51-S45
     19 1215  347 (S39W35 L=240)   360     53        1 (1182/ 322-1242/ 364)  14.0 W32-W37 S41-S37
  

    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