Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1017 1715 (N45E16 L=144)   291     39        0 ( 989/1701-1047/1729)  26.6 E17-E12 N43-N46
      2 1172 1628 (N41E00 L=159)   277     63        0 (1137/1603-1203/1647)  53.1 E02-W01 N38-N43
      3  403 1537 (N22E55 L=104)  1274    235        2 ( 371/1456- 427/1631) -63.4 E59-E51 N16-N28
      4  383 1405 (N13E51 L=108)   860    315        4 ( 309/1351- 501/1454)  29.7 E56-E42 N08-N17
      5  579 1302 (N10E33 L=127)   514    129        2 ( 544/1254- 618/1346) -59.0 E35-E29 N06-N12
      6 1120 1244 (N14W02 L=162)   233     49        0 (1087/1218-1142/1265) -26.6 W00-W04 N12-N15
      7  615  934 (S12E24 L=135)   723    196        3 ( 577/ 847- 668/1024) -63.4 E27-E19 S16-S06
      8 1308  947 (S01W18 L=178)   270     76        1 (1272/ 919-1352/ 974)  36.9 W16-W20 S02-N01
      9  245  628 (S38E57 L=102)   563    170        0 ( 222/ 600- 263/ 661)  39.8 E60-E53 S40-S35
     10 1317  404 (S35W34 L=193)   433    119        1 (1278/ 383-1374/ 434)   0.0 W30-W38 S35-S32
  

    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