Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  433 1539 (N13E55 L= 17)  1384    321        4 ( 371/1434- 524/1659) -82.9 E58-E51 N06-N22
      2 1312 1540 (N35W05 L= 76)  1396    291        2 (1231/1423-1359/1644) -49.8 E03-W10 N28-N41
      3  692 1446 (N13E31 L= 41)   254     79        0 ( 670/1407- 711/1480)  90.0 E31-E29 N10-N15
      4  380 1340 (N00E50 L= 22)   285     44        2 ( 354/1317- 406/1364)   0.0 E52-E47 N00-N01
      5 1622 1167 (N20W34 L=106)  1333    245        3 (1563/1061-1686/1261)  77.9 W30-W38 N13-N27
      6 1840 1041 (N19W57 L=129)  1569    167        3 (1806/ 966-1869/1134) -45.0 W52-W62 N14-N24
      7 1326 1025 (N04W17 L= 89)   362     80        0 (1281/1004-1358/1052)  36.9 W14-W19 N02-N06
      8  972  729 (S21W04 L= 76)   631    249        1 ( 872/ 689-1094/ 762)  39.3 E01-W10 S26-S16
      9  542  702 (S34E24 L= 47)  1019    200        1 ( 464/ 676- 647/ 747)   7.6 E31-E16 S35-S31
     10 1606  555 (S15W52 L=124)   534    109        1 (1571/ 516-1656/ 591)  56.3 W50-W55 S17-S11
     11 1674  549 (S13W59 L=130)   400     61        1 (1646/ 525-1697/ 578)  68.2 W57-W61 S15-S10
     12 1562  398 (S25W60 L=131)   523    114        1 (1503/ 377-1599/ 426)  39.8 W55-W64 S27-S22
  

    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