Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  593 1707 (N58E35 L= 63)  1723    318        2 ( 424/1679- 728/1739)  -1.2 E66-E17 N55-N60
      2  516 1364 (N33E33 L= 65)  1191    321        0 ( 428/1235- 582/1467) -50.5 E43-E26 N23-N40
      3 1799 1367 (N14W65 L=163)   533     61        2 (1772/1346-1830/1391) -26.6 W61-W68 N12-N15
      4  391 1273 (N28E43 L= 55)   405    119        1 ( 365/1231- 426/1330) -49.4 E46-E39 N24-N31
      5 1752 1215 (N06W53 L=151)   321     32        1 (1733/1186-1767/1239) -90.0 W51-W54 N04-N07
      6 1564  820 (S14W32 L=130)   351     90        0 (1542/ 773-1583/ 863)  59.0 W30-W33 S16-S11
      7 1432  539 (S31W22 L=120)  2536    529        4 (1196/ 405-1641/ 649)  15.9 W03-W38 S37-S26
      8  559  603 (S15E36 L= 62)   630    138        0 ( 495/ 571- 629/ 638)   9.5 E42-E30 S16-S12
      9  697  206 (S48E50 L= 48)  1021     76        2 ( 656/ 184- 743/ 228)   3.8 E59-E43 S50-S45
  

    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