Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  499 1682 (N62E30 L= 39)   514     36        1 ( 477/1657- 521/1704)   9.5 E34-E25 N59-N65
      2 1314 1573 (N28W35 L=103)   441    139        1 (1294/1512-1346/1632)  54.5 W32-W38 N23-N30
      3  889 1528 (N37W06 L= 74)   426    140        0 ( 825/1497- 950/1547) -40.6 W01-W09 N33-N39
      4  499 1487 (N46E26 L= 42)   871    194        0 ( 418/1443- 574/1513) -10.6 E35-E19 N42-N47
      5  813 1373 (N29E04 L= 65)   232     81        1 ( 782/1345- 844/1406)  11.3 E06-E01 N28-N31
      6  731 1265 (N25E12 L= 57)   241    145        0 ( 673/1245- 783/1296) -39.8 E14-E08 N23-N28
      7  461 1198 (N27E32 L= 37)   527    173        1 ( 431/1156- 523/1262) -53.1 E33-E27 N23-N31
      8  820 1141 (N16E09 L= 60)   214     91        1 ( 804/1090- 840/1182)  53.1 E10-E07 N12-N17
      9  293 1005 (N20E49 L= 19)   424    134        0 ( 258/ 958- 336/1065) -53.1 E51-E45 N15-N23
     10  658 1031 (N13E21 L= 47)   381     99        1 ( 629/ 984- 697/1061)  21.8 E23-E17 N10-N14
     11 1550 1022 (S09W31 L= 99)   237     76        0 (1525/ 984-1569/1058)  36.9 W28-W32 S10-S07
     12 1531  759 (S24W24 L= 93)   694    246        2 (1501/ 681-1584/ 882) -90.0 W22-W27 S30-S16
     13 1715  701 (S34W39 L=108)   499    101        2 (1693/ 652-1733/ 743)  80.5 W37-W41 S36-S30
     14 1305  643 (S25W07 L= 75)   435    149        1 (1219/ 626-1370/ 661) -21.8 W01-W11 S26-S22
     15  637  319 (S28E49 L= 19)  2217    310        2 ( 580/ 197- 751/ 429)  18.9 E72-E34 S36-S23
     16 1462  348 (S52W13 L= 81)   585     93        1 (1411/ 330-1504/ 370) -29.1 W08-W17 S55-S49
     17 1336  241 (S58E03 L= 65)   749     73        0 (1305/ 217-1370/ 263)  20.6 E07-W01 S60-S55
  

    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