Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1376 1429 (N18W33 L= 54)   916    278        2 (1261/1407-1521/1450) -31.0 W26-W41 N12-N21
      2 1731 1347 (N02W55 L= 76)   981    284        2 (1674/1286-1802/1403)   6.3 W47-W65 N00-N03
      3 1422 1086 (S02W24 L= 45)   399    114        0 (1402/1038-1456/1145)  71.6 W22-W26 S05-N01
      4 1540 1001 (S10W29 L= 50)   615    211        1 (1455/ 972-1632/1050) -17.1 W23-W36 S12-S06
      5 1768  954 (S20W45 L= 66)   338     72        1 (1742/ 925-1797/ 975)   0.0 W42-W47 S20-S18
      6 1595  818 (S22W28 L= 49)   915    220        1 (1562/ 731-1628/ 951)  59.7 W25-W33 S27-S15
      7 1885  866 (S30W55 L= 76)  1614    297        0 (1825/ 818-1944/ 941)   8.1 W48-W64 S32-S27
      8  644  791 (N03E28 L=353)   449    116        1 ( 607/ 734- 679/ 841)  29.7 E31-E24 N00-N05
      9  535  745 (N03E36 L=345)   331     88        0 ( 493/ 726- 580/ 772) -38.7 E38-E33 N00-N04
     10 1263  721 (S17W05 L= 26)   220     47        1 (1244/ 692-1282/ 751) -90.0 W04-W06 S19-S15
  

    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