Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  808 1625 (N23E32 L=206)   445    142        1 ( 778/1561- 851/1672) -83.7 E33-E30 N18-N27
      2  526 1615 (N16E50 L=188)   345    108        1 ( 501/1584- 552/1644) -26.6 E53-E47 N14-N17
      3  456 1482 (N06E49 L=189)   347     72        0 ( 425/1465- 494/1496)   9.5 E52-E46 N05-N07
      4 1368 1347 (N22W11 L=249)   417     60        1 (1338/1314-1396/1382) -36.9 W09-W13 N20-N24
      5  302 1170 (S15E51 L=188)  1336    330        2 ( 202/1083- 376/1298)  59.0 E58-E44 S21-S06
      6  202 1208 (S15E62 L=176)  1063    167        2 ( 166/1152- 249/1256)  81.9 E65-E57 S18-S11
      7  279 1198 (S14E54 L=185)   461    112        1 ( 248/1163- 307/1223)  63.4 E56-E51 S16-S12
      8  664 1048 (S14E21 L=217)   446    201        0 ( 581/ 965- 708/1113) -29.1 E27-E17 S17-S10
      9 1042  852 (S14W06 L=244)   458    114        2 ( 990/ 789-1096/ 903)  71.6 W04-W07 S19-S10
     10  197  790 (S40E56 L=182)   663    161        0 ( 165/ 770- 237/ 822)  -8.1 E59-E51 S41-S38
     11 1010  666 (S26W10 L=248)   251     64        0 ( 984/ 637-1039/ 704) -21.8 W07-W12 S26-S24
     12 1431  667 (S13W35 L=273)   400    100        0 (1389/ 623-1482/ 703)  66.8 W32-W36 S16-S09
     13  831  651 (S32E01 L=238)   258     53        1 ( 798/ 625- 859/ 669) -14.0 E03-W01 S32-S30
     14 1709  477 (S13W67 L=305)   629    122        2 (1688/ 440-1725/ 519) -45.0 W63-W70 S15-S11
     15  848  388 (S49W11 L=250)   303     90        0 ( 831/ 354- 861/ 427) -71.6 W09-W13 S51-S45
     16  972  295 (S51W27 L=265)  1196    260        6 ( 926/ 219-1016/ 379) -38.7 W20-W35 S56-S44
     17 1273  333 (S37W43 L=281)   697    100        2 (1235/ 305-1317/ 366)  45.0 W38-W45 S39-S34
     18 1117  282 (S46W37 L=275)   418    104        2 (1074/ 253-1161/ 297)  31.0 W33-W40 S47-S43
  

    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