Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1041 1820 (N54W13 L=347)   461     34        0 (1015/1798-1067/1844)  26.6 W10-W16 N51-N56
      2  623 1702 (N48E28 L=306)   450     51        1 ( 606/1670- 643/1731) -56.3 E30-E25 N45-N51
      3 1371 1564 (N29W30 L=  4)   451     96        3 (1348/1510-1401/1601)  45.0 W27-W33 N25-N31
      4  399 1446 (N31E43 L=291)   422     66        1 ( 371/1420- 430/1472)  21.8 E45-E40 N29-N32
      5  490 1407 (N27E34 L=300)   385     72        0 ( 451/1393- 535/1423)   0.0 E37-E30 N26-N28
      6  433  603 (S21E46 L=288)   361     69        1 ( 407/ 574- 468/ 630)  18.4 E48-E42 S22-S19
      7  709  463 (S33E29 L=305)   711    192        2 ( 643/ 419- 797/ 519)  17.4 E36-E20 S35-S29
      8  390  389 (S35E66 L=268)   666    139        1 ( 339/ 363- 426/ 429)   6.0 E75-E56 S36-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