Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  834 1819 (N62W19 L=113)   366     70        0 ( 807/1795- 866/1842)  10.3 W14-W25 N60-N63
      2 1159 1806 (N47W43 L=137)   347     41        1 (1140/1783-1184/1826)   0.0 W39-W46 N45-N48
      3 1411 1742 (N34W55 L=149)   932     93        3 (1364/1717-1458/1766) -63.4 W50-W58 N31-N37
      4  968 1708 (N47W20 L=115)   257     52        0 ( 939/1687- 996/1722) -53.1 W18-W21 N44-N48
      5  801 1615 (N46W03 L= 97)   244     63        1 ( 782/1583- 826/1642) -80.5 W01-W04 N42-N48
      6  997 1383 (N25W08 L=102)   195     84        0 ( 975/1342-1022/1423) -90.0 W07-W08 N22-N28
      7  563 1400 (N38E21 L= 73)   296     81        1 ( 526/1376- 601/1419) -26.6 E24-E18 N35-N39
      8  272 1393 (N45E53 L= 41)   688     75        1 ( 231/1376- 312/1415) -13.0 E60-E47 N43-N46
      9  513 1243 (N29E27 L= 67)   220     99        0 ( 496/1211- 544/1267)  45.0 E28-E24 N27-N30
     10  223 1105 (N27E56 L= 39)  2097    258        7 ( 187/1019- 284/1218) -58.4 E60-E49 N21-N34
     11  397 1122 (N25E38 L= 56)   429    106        1 ( 357/1100- 449/1150)   0.0 E41-E33 N23-N26
     12  285  881 (N13E52 L= 42)   578     97        3 ( 256/ 837- 307/ 923)  68.2 E54-E49 N10-N15
     13 1178  852 (S09W04 L= 99)   408    155        0 (1153/ 774-1213/ 916)  45.0 W01-W07 S12-S06
     14 1758  655 (S39W45 L=139)   989    150        1 (1695/ 624-1828/ 679)  -3.6 W37-W54 S40-S35
     15  934  484 (S24E20 L= 74)   297    119        1 ( 875/ 471- 991/ 502) -14.0 E24-E16 S25-S22
     16  824  407 (S27E31 L= 64)   246    100        1 ( 799/ 365- 856/ 440)  39.8 E33-E27 S29-S24
     17  477  300 (S26E68 L= 26)   735     75        0 ( 438/ 282- 513/ 320)   5.7 E81-E61 S27-S25
     18 1253  254 (S52E10 L= 84)   801    109        1 (1207/ 233-1312/ 276) -24.0 E14-E04 S54-S49
     19 1140  175 (S56E27 L= 68)   542     68        0 (1112/ 156-1167/ 200)  20.6 E30-E21 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