Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  450 1432 (N33E35 L=189)   645     93        3 ( 401/1400- 495/1469)  23.2 E39-E31 N31-N35
      2  295 1418 (N35E49 L=175)   364     71        1 ( 264/1395- 317/1437)  11.3 E52-E46 N33-N35
      3  645 1265 (N19E20 L=204)   416    123        0 ( 607/1197- 691/1306)  45.0 E23-E17 N15-N21
      4 1203 1269 (N12W14 L=238)   221     60        0 (1172/1228-1229/1296) -63.4 W12-W15 N09-N13
      5  560  887 (S02E30 L=194)   231     72        0 ( 529/ 846- 586/ 920) -63.4 E31-E28 S04-N00
      6  297  706 (S09E55 L=169)  1232    312        4 ( 255/ 598- 333/ 819)  55.3 E63-E50 S14-S01
      7  478  562 (S20E44 L=180)   724    223        2 ( 430/ 489- 524/ 644)  39.3 E50-E37 S24-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