Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1775 1479 (N19W71 L= 10)   727     87        2 (1747/1447-1793/1519)  12.1 W65-W79 N16-N20
      2 1377 1448 (N26W32 L=331)   291     63        1 (1350/1429-1410/1469) -45.0 W29-W33 N24-N27
      3 1193 1357 (N24W16 L=315)   292    153        0 (1164/1292-1213/1421)  77.5 W14-W17 N19-N28
      4 1794 1391 (N13W67 L=  6)   588     41        1 (1773/1360-1811/1420)  45.0 W63-W69 N11-N15
      5 1538 1321 (N15W39 L=338)   398    111        0 (1478/1289-1582/1348) -39.8 W35-W42 N12-N17
      6  726 1232 (N23E16 L=283)   419     97        3 ( 691/1204- 773/1258) -21.8 E18-E13 N21-N24
      7  269 1046 (N16E54 L=245)   451    123        0 ( 255/ 993- 290/1098)  74.1 E55-E52 N12-N19
      8  522  757 (S04E36 L=263)   360    127        0 ( 457/ 739- 572/ 772) -12.5 E41-E32 S05-S02
      9  755  398 (S30E31 L=269)  2404    414        5 ( 605/ 295- 882/ 521)  22.1 E49-E17 S37-S22
     10 1315  428 (S38W11 L=310)   433     95        0 (1280/ 394-1350/ 460)  38.7 W08-W14 S39-S35
  

    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