Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1038 1770 (N40E22 L=263)   298     46        0 (1008/1758-1068/1784)  45.0 E24-E20 N38-N41
      2 1121 1713 (N38E14 L=271)   676    109        2 (1072/1682-1173/1748)   8.1 E16-E09 N35-N40
      3  768 1641 (N23E33 L=252)   481    168        1 ( 743/1563- 793/1709) -65.6 E35-E30 N17-N28
      4  544 1522 (N11E43 L=242)  2089    486        3 ( 481/1377- 662/1673) -32.0 E54-E30 N04-N20
      5 1078 1328 (N11E04 L=281)  1225    382        2 (1041/1158-1110/1530) -70.0 E09-W01 N01-N23
      6  112 1122 (S18E71 L=214)  4797    404        1 (  74/1046- 146/1205)  83.7 E76-E66 S22-S13
      7 1814 1044 (N16W51 L=336)  1135    160        1 (1773/ 973-1848/1111) -31.0 W45-W55 N12-N18
      8  312  982 (S24E46 L=239)  1139    264        0 ( 227/ 913- 398/1052)  37.6 E53-E39 S29-S19
      9 1813  832 (N04W57 L=342)   701    145        2 (1774/ 796-1850/ 870)  35.5 W53-W60 N01-N07
     10 1433  595 (S20W37 L=322)   306     74        0 (1408/ 555-1451/ 637) -78.7 W35-W38 S22-S17
  

    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