Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  757 1810 (N42E40 L= 92)  1421    138        3 ( 707/1755- 805/1852) -17.4 E49-E33 N37-N45
      2 1023 1684 (N36E12 L=120)   378     73        2 ( 995/1649-1052/1717) -31.0 E14-E09 N33-N38
      3  493 1647 (N26E51 L= 81)   650    118        0 ( 446/1605- 543/1673) -15.9 E57-E43 N23-N27
      4 1538 1616 (N42W30 L=162)   526     77        2 (1514/1586-1572/1650)  40.6 W26-W33 N38-N44
      5 1721 1243 (N19W43 L=175)   509    112        0 (1698/1211-1750/1291) -51.3 W41-W46 N17-N22
      6  132  992 (S17E70 L= 62)  4503    248       11 (  86/ 891- 210/1078)  21.4 E84-E58 S22-S12
      7  633 1059 (S09E24 L=108)   214     57        0 ( 603/1040- 668/1076)  26.6 E26-E21 S10-S08
      8 1871  940 (N05W62 L=194)   501     78        0 (1847/ 899-1888/ 977) -33.7 W58-W65 N02-N06
      9 1512  909 (S05W32 L=164)   271     92        0 (1458/ 894-1554/ 921)   8.1 W27-W34 S05-S03
     10  418  881 (S22E39 L= 93)   293     65        3 ( 394/ 852- 443/ 908) -56.3 E40-E36 S24-S20
     11  900  874 (S17E05 L=127)   299     78        0 ( 857/ 857- 932/ 904)  45.0 E07-E03 S18-S14
     12  174  782 (S30E68 L= 64)  1254    171        2 ( 155/ 701- 193/ 869)  79.7 E71-E64 S35-S24
     13 1532  678 (S18W39 L=171)   379     72        0 (1497/ 645-1564/ 713)  63.4 W37-W40 S19-S15
     14 1570  571 (S23W47 L=178)   430     81        0 (1528/ 544-1607/ 596)  -7.1 W42-W50 S23-S22
     15  607  424 (S51E26 L=106)   788    107        1 ( 581/ 364- 631/ 473) -77.5 E28-E23 S56-S47
     16 1533  438 (S32W51 L=183)   335     63        1 (1508/ 409-1557/ 460) -45.0 W48-W53 S33-S30
  

    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