Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1339 1673 (N46W32 L=212)   949    260        1 (1239/1617-1444/1736) -38.2 W25-W39 N40-N52
      2 1666 1501 (N31W56 L=235)   555     94        1 (1620/1478-1715/1521) -18.4 W50-W60 N28-N32
      3 1231 1487 (N33W17 L=197)   376     82        0 (1193/1462-1276/1511) -31.0 W14-W20 N30-N34
      4  547 1294 (N22E32 L=148)   567    151        1 ( 494/1237- 583/1322) -15.9 E36-E29 N18-N23
      5 1249 1251 (N17W15 L=195)   413    160        0 (1195/1193-1290/1318) -58.0 W12-W17 N13-N21
      6  384 1216 (N17E45 L=135)   272     32        0 ( 370/1192- 400/1240)  90.0 E46-E43 N15-N18
      7  879 1056 (N07E09 L=171)   262     83        0 ( 839/1042- 928/1068)   0.0 E11-E05 N05-N07
      8 1101  865 (S06W04 L=184)   709    190        2 (1020/ 805-1165/ 923) -45.0 W00-W08 S09-S01
      9 1541  806 (S12W33 L=213)   932    287        0 (1493/ 722-1584/ 909) -80.5 W29-W36 S16-S04
     10 1611  806 (S12W39 L=219)   264     69        1 (1587/ 772-1625/ 829) -63.4 W36-W40 S14-S10
     11  253  744 (S13E59 L=121)  1120    182        3 ( 232/ 652- 285/ 817)  57.5 E64-E54 S18-S07
     12 1154  612 (S22W07 L=187)  1431    360        2 (1046/ 495-1281/ 736)  46.8 W00-W15 S30-S14
     13  756  685 (S16E18 L=161)  1042    242        4 ( 649/ 642- 850/ 726)  -8.7 E25-E12 S19-S13
     14 1546  602 (S25W36 L=216)   574    164        0 (1497/ 538-1598/ 663)  53.1 W32-W39 S29-S21
     15  847  586 (S23E13 L=166)   587    121        4 ( 805/ 534- 886/ 624)  45.0 E16-E10 S26-S20
     16 1002  231 (S54E07 L=173)   548     54        1 ( 978/ 203-1028/ 260) -63.4 E08-E04 S56-S50
     17 1120  204 (S58W06 L=186)  1075    128        2 (1064/ 177-1168/ 233)   0.0 W00-W12 S61-S53
     18 1222  188 (S61W20 L=199)   738     54        0 (1204/ 158-1239/ 224) -78.7 W16-W24 S65-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