Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  988 1754 (N48W12 L=164)   348     32        0 ( 968/1729-1011/1775)  21.8 W09-W14 N46-N50
      2  869 1574 (N36E02 L=150)   226     38        0 ( 841/1559- 899/1587) -14.0 E04-W00 N35-N37
      3  760 1549 (N36E10 L=142)   358     92        0 ( 727/1509- 793/1585)  38.7 E13-E07 N34-N38
      4 1389  885 (S13W20 L=173)   215     37        0 (1371/ 854-1408/ 912)  56.3 W18-W21 S14-S11
      5 1311  727 (S22W14 L=166)  1195    378        0 (1148/ 668-1450/ 820)   5.7 W03-W23 S24-S17
      6  421  647 (S14E46 L=106)   707     96        1 ( 381/ 600- 451/ 693)  45.0 E49-E43 S16-S11
      7  853  596 (S23E17 L=135)   493    170        0 ( 825/ 509- 891/ 665)  48.4 E21-E13 S28-S19
      8  486  637 (S15E41 L=111)   259     43        0 ( 466/ 609- 501/ 659) -90.0 E42-E40 S17-S13
      9  223  559 (S16E74 L= 78)  2314    137        1 ( 183/ 497- 268/ 612)  69.4 E81-E67 S20-S12
     10  547  389 (S32E48 L=104)  1755    290        0 ( 468/ 336- 630/ 472) -45.0 E55-E39 S35-S25
  

    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