Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  988 1727 (N38E16 L=298)   523     88        1 ( 935/1696-1026/1751)  -8.1 E20-E13 N35-N39
      2  727 1631 (N26E31 L=284)   287     58        1 ( 707/1605- 755/1654)  90.0 E32-E28 N24-N28
      3  235 1379 (N06E63 L=252)  1352    142        2 ( 177/1329- 279/1431)  36.9 E68-E58 N02-N09
      4 1344 1342 (N17W14 L=329)   259     72        0 (1302/1325-1387/1356)   0.0 W11-W17 N16-N18
      5  541 1284 (N02E33 L=281)   608    126        2 ( 507/1224- 582/1329)  63.4 E36-E31 S01-N05
      6 1716 1209 (N17W42 L=357)  1653    235        5 (1671/1124-1785/1312) -70.0 W38-W48 N11-N22
      7  390 1278 (N00E44 L=270)   276     67        0 ( 352/1265- 426/1291)   0.0 E47-E41 N00-N01
      8 1640 1158 (N12W36 L=351)   521    118        3 (1605/1112-1668/1212)  90.0 W34-W38 N09-N15
      9 1391 1171 (N08W19 L=334)   271     61        1 (1370/1135-1415/1202)  78.7 W18-W20 N05-N10
     10 1488 1141 (N08W26 L=340)   231     41        1 (1466/1115-1512/1165) -18.4 W24-W27 N06-N09
     11 1860  831 (S02W62 L= 17)   419     52        0 (1826/ 817-1889/ 847)  14.0 W58-W66 S02-N00
     12  649  780 (S26E20 L=294)  1375    415        2 ( 501/ 715- 860/ 826)  17.7 E31-E06 S31-S21
     13 1462  586 (S24W37 L=352)  3687    830        6 (1185/ 509-1755/ 721)   3.4 W15-W66 S29-S13
     14  474  524 (S45E37 L=278)   859    175        3 ( 378/ 498- 559/ 551)   0.0 E47-E28 S46-S43
     15  704  494 (S44E15 L=300)   398     99        1 ( 658/ 478- 754/ 518)  -7.1 E18-E10 S45-S42
     16  875  468 (S43E00 L=314)   295     91        0 ( 831/ 450- 921/ 492)   0.0 E04-W03 S43-S41
     17 1098  407 (S43W19 L=333)   312    111        0 (1053/ 387-1152/ 439)  -6.3 W14-W23 S43-S41
  

    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