Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1070 1868 (N54E31 L= 28)   525     57        1 (1044/1848-1094/1891)  -8.1 E35-E27 N51-N56
      2  529 1770 (N30E61 L=358)   631    165        0 ( 514/1748- 553/1798) -20.6 E65-E57 N28-N31
      3  938 1576 (N27E21 L= 38)   863    248        1 ( 850/1495-1048/1657) -12.5 E29-E11 N23-N30
      4 1342 1556 (N37W06 L= 65)   421     68        2 (1302/1528-1382/1583)   8.1 W02-W09 N35-N38
      5  324 1494 (N08E62 L=357)  1371    220        3 ( 262/1428- 381/1552)  48.8 E68-E58 N02-N12
      6  485 1299 (S01E40 L= 19)   570    220        2 ( 412/1277- 598/1325)  17.1 E45-E32 S02-N03
      7  287 1175 (S12E53 L=  6)   333     72        0 ( 270/1139- 301/1200) -45.0 E54-E51 S13-S10
      8 1763  773 (N02W56 L=115)   692    144        0 (1737/ 733-1801/ 811) -16.7 W51-W61 N00-N03
  

    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