Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  661 1569 (N48E19 L=224)   364     75        2 ( 623/1548- 696/1588) -26.6 E23-E15 N46-N50
      2  817 1554 (N44E05 L=238)   317     59        1 ( 780/1532- 846/1573) -21.8 E07-E02 N41-N45
      3 1343 1372 (N21W27 L=270)  1499    534        0 (1105/1312-1547/1434) -30.7 W12-W39 N13-N29
      4  963 1373 (N29W02 L=245)   279    144        0 ( 900/1337-1008/1392) -31.0 E01-W04 N25-N30
      5  813 1290 (N26E09 L=234)   329    103        2 ( 774/1258- 857/1328) -45.0 E11-E06 N23-N29
      6  314 1259 (N31E51 L=192)   336    101        0 ( 276/1224- 358/1280)  14.0 E54-E45 N28-N32
      7 1298 1222 (N13W20 L=263)   404    179        1 (1212/1193-1373/1246) -33.7 W15-W24 N10-N16
      8  344 1154 (N24E47 L=196)  1532    226        4 ( 309/1056- 386/1241)  80.5 E50-E42 N17-N29
      9 1443 1100 (N03W27 L=270)   275     64        0 (1425/1065-1457/1132)  63.4 W25-W28 N01-N05
     10 1853 1095 (S06W62 L=305)   506     43        0 (1833/1071-1874/1114)  21.8 W58-W65 S07-S05
     11 1651  908 (S12W39 L=282)   253     62        0 (1620/ 887-1677/ 925) -26.6 W36-W40 S14-S11
     12  226  866 (N07E61 L=182)   615     67        1 ( 205/ 839- 249/ 896)  36.9 E64-E57 N05-N08
     13  294  648 (S07E60 L=183)   632    131        1 ( 250/ 602- 349/ 693) -45.0 E64-E56 S10-S03
     14 1104  531 (S25E04 L=239)   363    121        0 (1056/ 494-1158/ 562)  23.2 E07-W00 S27-S23
     15  559  402 (S26E47 L=196)   949    203        1 ( 519/ 320- 627/ 462)  24.2 E58-E38 S32-S22
     16  708  330 (S34E39 L=204)  1027    203        2 ( 641/ 279- 766/ 398)  21.3 E49-E31 S37-S28
     17 1406  289 (S54W18 L=261)   414     46        0 (1378/ 273-1436/ 301) -15.9 W14-W21 S56-S51
  

    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