Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1401 1653 (N39W27 L=256)  2173    337        2 (1271/1581-1536/1707) -18.4 W16-W37 N34-N43
      2  887 1535 (N27E11 L=217)   301     67        0 ( 850/1514- 926/1555) -11.3 E13-E08 N26-N28
      3  316 1506 (N25E58 L=171)  1142    102        3 ( 277/1457- 354/1552) -49.4 E64-E53 N21-N28
      4 1669 1418 (N24W45 L=274)   576     84        0 (1650/1378-1698/1466) -80.5 W42-W47 N21-N27
      5 1495 1178 (N07W29 L=258)   344     83        1 (1475/1143-1512/1225)  90.0 W27-W30 N05-N10
      6  484  936 (S11E34 L=194)   449     97        0 ( 448/ 891- 520/ 988)  56.3 E37-E31 S13-S07
      7  254  877 (S14E55 L=174)  1621    242        4 ( 216/ 757- 293/ 979)  85.9 E58-E51 S21-S07
      8  901  783 (S19E07 L=222)   311    127        1 ( 860/ 756- 960/ 837)  39.8 E09-E03 S20-S15
      9 1703  761 (S16W49 L=278)   474     84        3 (1678/ 723-1736/ 798) -33.7 W45-W52 S18-S13
     10  449  623 (S30E41 L=187)  1191    186        3 ( 415/ 540- 491/ 710) -84.8 E43-E38 S35-S24
     11 1719  639 (S23W55 L=283)   486     68        1 (1692/ 610-1743/ 672)  45.0 W51-W56 S25-S21
     12 1352  600 (S29W25 L=254)   688    119        3 (1308/ 561-1386/ 653)  63.4 W21-W27 S31-S25
     13 1105  591 (S31W07 L=236)   640    206        0 (1009/ 572-1203/ 621)   8.7 W00-W13 S31-S28
     14  647  535 (S36E26 L=202)   416    108        1 ( 589/ 521- 704/ 549)   0.0 E31-E21 S36-S34
     15  848  501 (S38E11 L=218)   290     54        2 ( 819/ 470- 876/ 532)  59.0 E13-E08 S40-S35
     16  568  280 (S56E50 L=178)   564     52        1 ( 540/ 258- 598/ 307)  11.3 E59-E44 S58-S53
  

    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