Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  838 1808 (N38E43 L=217)  1358    193        1 ( 759/1780- 933/1841)   3.2 E53-E34 N36-N40
      2  484 1690 (N21E62 L=199)   498     60        0 ( 468/1668- 505/1718) -63.4 E65-E57 N19-N23
      3  730 1641 (N24E38 L=223)   304     50        2 ( 703/1618- 760/1669)  90.0 E39-E36 N21-N25
      4  296 1538 (N08E70 L=190)  1490    148        3 ( 249/1472- 344/1581)  66.0 E75-E65 N02-N11
      5 1288 1369 (N22W06 L=266)   484    193        1 (1247/1328-1346/1427) -14.0 W02-W10 N20-N25
      6  202 1296 (S08E68 L=193)   619     62        1 ( 169/1277- 234/1317)  18.4 E72-E63 S09-S06
      7 1819  855 (N09W58 L=319)  1418    229        2 (1763/ 764-1860/ 950)  77.0 W55-W63 N02-N15
      8  722  572 (S40E06 L=254)   383    116        0 ( 700/ 522- 759/ 619) -39.8 E08-E02 S41-S36
      9 1357  535 (S22W36 L=296)   336    151        0 (1311/ 515-1424/ 572)  14.0 W31-W39 S23-S19
  

    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