Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  370 1085 (N24E42 L=317)   627    159        1 ( 308/1049- 434/1117) -26.6 E47-E37 N21-N26
      2  306  968 (N18E49 L=310)   671    146        1 ( 255/ 929- 357/1013) -36.9 E53-E45 N15-N21
      3 1371  910 (S08W17 L= 17)   419     78        1 (1322/ 877-1417/ 935)   0.0 W13-W20 S08-S06
      4 1477  898 (S11W24 L= 23)   225     44        0 (1461/ 871-1496/ 929)  71.6 W22-W24 S12-S09
      5  654  827 (N04E26 L=333)   515     83        0 ( 625/ 789- 687/ 876) -80.5 E27-E24 N01-N07
      6 1193  761 (S12W04 L=  3)   764    195        0 (1149/ 667-1223/ 858)  55.0 E01-W06 S16-S06
      7 1655  785 (S23W35 L= 34)   261     72        0 (1636/ 748-1678/ 826)  53.1 W32-W36 S25-S21
      8 1308  762 (S15W10 L= 10)   220     36        0 (1284/ 743-1329/ 784)  14.0 W08-W12 S15-S13
      9  652  728 (S02E29 L=330)   258     46        0 ( 629/ 708- 679/ 756) -45.0 E30-E27 S03-N00
     10 1185  517 (S27E03 L=356)   796    245        0 (1075/ 432-1293/ 611)  18.4 E12-W06 S30-S23
     11  837  432 (S24E27 L=332)   312     68        0 ( 810/ 402- 866/ 464)  31.0 E29-E24 S25-S22
     12  832  236 (S39E39 L=320)   797    238        1 ( 808/ 174- 858/ 298)  45.0 E44-E32 S45-S34
  

    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