Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  752 1656 (N37E22 L=257)  1508    289        2 ( 638/1581- 889/1728) -24.6 E34-E10 N31-N43
      2  957 1416 (N20E05 L=275)   209     51        1 ( 935/1393- 975/1448) -76.0 E06-E03 N18-N22
      3 1904 1093 (N04W66 L=346)   595     77        0 (1884/1053-1921/1121) -90.0 W63-W69 N01-N05
      4 1838 1048 (N01W58 L=338)   398     78        1 (1818/1009-1860/1079)  45.0 W55-W60 S01-N02
      5 1719  885 (S10W47 L=327)   351     67        1 (1699/ 855-1742/ 926) -36.9 W45-W49 S11-S07
      6 1541  819 (S15W34 L=314)   452    133        0 (1503/ 759-1559/ 874) -60.3 W30-W35 S18-S11
      7  304  773 (S19E52 L=228)  1314    239        3 ( 268/ 657- 343/ 847)  50.7 E58-E47 S25-S14
      8 1710  686 (S22W51 L=331)   955    220        3 (1657/ 599-1737/ 782)  90.0 W48-W55 S27-S15
      9  627  595 (S31E28 L=252)   639    137        1 ( 585/ 541- 661/ 658)  58.0 E31-E24 S34-S26
     10  423  575 (S32E46 L=234)   587    113        2 ( 379/ 549- 481/ 599)  14.0 E51-E39 S33-S30
     11 1333  217 (S59W42 L=322)   425     42        1 (1308/ 200-1356/ 236)  26.6 W37-W47 S60-S57
     12 1174  188 (S63W24 L=304)   794     94        1 (1130/ 158-1205/ 224)  60.3 W19-W28 S67-S59
  

    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