Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  593 1733 (N41E43 L=295)   527     76        0 ( 552/1712- 633/1751) -11.3 E47-E37 N39-N42
      2  498 1706 (N39E51 L=287)   592     63        0 ( 477/1669- 524/1736)  90.0 E53-E48 N35-N41
      3 1574 1563 (N34W39 L= 17)   309     31        0 (1557/1535-1587/1584) -90.0 W37-W40 N31-N35
      4  341 1545 (N27E57 L=280)   608     64        1 ( 319/1515- 367/1572)  90.0 E60-E53 N24-N28
      5  334  841 (S17E47 L=290)   386     96        0 ( 321/ 786- 350/ 889)  63.4 E48-E45 S20-S14
      6  221  762 (S22E61 L=276)   994     86        1 ( 169/ 737- 273/ 791)   3.8 E69-E54 S23-S20
      7 1688  522 (S30W58 L= 35)  1430    159        4 (1636/ 478-1762/ 563)  29.1 W52-W66 S32-S26
      8  532  508 (S39E37 L=301)  1031    116        2 ( 479/ 474- 585/ 543)  14.0 E43-E31 S41-S36
      9 1194  444 (S40W17 L=355)   417    197        0 (1106/ 416-1272/ 495)  28.6 W10-W22 S42-S35
     10 1506  447 (S37W43 L= 21)   360     81        0 (1484/ 402-1522/ 484) -71.6 W41-W45 S40-S34
  

    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