Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1122 1873 (N58E30 L=219)   892     81        1 (1099/1831-1143/1918) -42.9 E39-E24 N52-N65
      2  765 1735 (N34E42 L=207)   370     74        0 ( 750/1698- 785/1772) -39.8 E45-E37 N31-N36
      3 1294 1652 (N43W00 L=249)   515     48        0 (1270/1622-1315/1685) -45.0 E02-W02 N41-N45
      4 1599 1582 (N49W30 L=279)   781    108        3 (1550/1552-1641/1613)  -5.2 W24-W36 N46-N51
      5 1684 1497 (N45W40 L=288)   385     55        1 (1666/1467-1703/1529) -78.7 W37-W42 N42-N47
      6 1637 1252 (N26W34 L=283)  1922    496        3 (1494/1057-1774/1415) -25.8 W20-W49 N17-N32
      7  199 1026 (S22E61 L=188)   834     73        1 ( 172/ 988- 228/1063)  51.3 E64-E56 S24-S19
      8 1084  591 (S27W16 L=264)  1428    399        5 ( 913/ 513-1252/ 641)  41.6 W06-W24 S35-S19
      9 1333  520 (S24W34 L=283)   312    129        0 (1280/ 476-1371/ 561)  82.9 W32-W35 S28-S20
     10 1440  340 (S32W55 L=303)  2121    326        2 (1263/ 316-1557/ 372)  23.5 W41-W65 S38-S27
     11  982  160 (S62W44 L=293)  1546    153        3 ( 948/ 103-1030/ 220)  -8.7 W31-W72 S66-S56
  

    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