Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  956 1592 (N40E07 L=342)   452     66        1 ( 929/1552- 983/1622)  90.0 E08-E04 N37-N42
      2 1691 1383 (N26W50 L= 38)  3282    764       12 (1464/1263-1897/1460) -15.4 W31-W71 N17-N31
      3  331 1413 (N25E54 L=294)   392    123        0 ( 306/1379- 360/1450) -29.7 E57-E50 N23-N27
      4  978 1328 (N22E04 L=345)   531    159        1 ( 941/1255-1011/1405) -68.2 E06-E01 N17-N27
      5 1700 1200 (N14W47 L= 36)   997    240        1 (1605/1111-1767/1290) -42.5 W40-W53 N08-N19
      6 1317 1082 (N07W18 L=  7)   456    149        1 (1292/1009-1355/1156) -71.6 W16-W20 N02-N11
      7 1463  985 (N01W28 L= 16)   808    197        3 (1376/ 946-1529/1037)  26.6 W22-W32 S01-N04
      8 1217  908 (S04W12 L=  0)   309    147        0 (1135/ 887-1277/ 927)   0.0 W07-W15 S04-S02
  

    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