Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1632 1645 (N22W64 L=102)  2781    314        2 (1521/1552-1725/1777)  74.5 W58-W71 N14-N32
      2 1322 1629 (N30W37 L= 75)  1353    243        2 (1217/1576-1430/1664)  -6.0 W27-W46 N28-N31
      3 1569 1571 (N20W52 L= 90)  2564    348        3 (1409/1458-1705/1631) -43.2 W42-W59 N11-N27
      4 1338 1485 (N21W31 L= 69)   323    105        0 (1285/1470-1390/1498) -23.2 W27-W34 N19-N22
      5 1711 1107 (S09W44 L= 82)  1383    287        2 (1645/1026-1755/1220)  70.0 W37-W48 S12-S01
      6  412 1119 (N22E36 L=  1)   639    136        1 ( 349/1067- 453/1153) -41.2 E41-E33 N18-N25
      7 1390  953 (S09W19 L= 57)   285     96        1 (1369/ 905-1411/1001)  80.5 W17-W20 S11-S05
      8  863  857 (S02E13 L= 25)  1589    385        6 ( 745/ 686- 952/ 982) -86.8 E17-E10 S13-N05
      9 1255  896 (S09W10 L= 48)  1212    260        7 (1134/ 826-1343/ 945)   7.6 W01-W16 S11-S06
     10  341  891 (N10E46 L=352)   616    128        1 ( 289/ 862- 405/ 930) -32.0 E50-E41 N07-N13
     11 1408  666 (S26W15 L= 53)  1515    359        3 (1298/ 618-1531/ 776) -40.2 W09-W24 S29-S17
     12  541  684 (S05E37 L=  1)   247     76        0 ( 523/ 659- 565/ 728) -90.0 E38-E35 S06-S02
     13 1116  397 (S36E11 L= 26)  3922    674       15 ( 958/ 184-1219/ 673)  42.3 E35-W00 S50-S18
     14  577  448 (S20E44 L=354)  1288    204        2 ( 523/ 376- 627/ 560)  35.5 E52-E37 S23-S13
     15  442  522 (S13E51 L=347)   401     87        3 ( 418/ 484- 467/ 558)  53.1 E54-E48 S14-S10
     16  437  405 (S20E59 L=339)   942     96        1 ( 406/ 363- 479/ 437)  10.3 E65-E54 S21-S17
  

    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