Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  666 1785 (N52E50 L=148)   604    104        0 ( 634/1741- 691/1812) -20.6 E58-E42 N48-N54
      2 1202 1472 (N32W10 L=208)  1587    463        2 (1070/1332-1374/1566) -29.7 W00-W21 N23-N38
      3  500 1388 (N21E40 L=158)   628    153        0 ( 474/1326- 547/1459) -41.6 E44-E34 N17-N25
      4 1688 1321 (N24W48 L=247)   366     61        0 (1673/1284-1704/1352) -90.0 W47-W50 N21-N26
      5 1686 1158 (N14W46 L=244)   428     84        1 (1654/1129-1721/1184) -18.4 W42-W48 N12-N15
      6 1858 1117 (N12W64 L=263)   901    115        1 (1835/1062-1884/1174)  49.4 W60-W69 N08-N15
      7  870 1094 (N06E10 L=188)   531    219        0 ( 819/1054- 919/1147) -45.0 E13-E07 N03-N09
      8 1750  859 (S04W53 L=251)   567     44        0 (1723/ 837-1773/ 885)  45.0 W50-W55 S05-S02
      9 1107  731 (S16W07 L=205)   411    139        0 (1061/ 691-1131/ 785)  81.9 W04-W08 S18-S11
     10 1294  643 (S20W21 L=219)   274     84        0 (1266/ 611-1326/ 685)  68.2 W18-W22 S22-S17
     11 1823  654 (S17W70 L=268)   672     54        0 (1793/ 633-1844/ 678)  -5.2 W64-W75 S18-S15
     12 1675  609 (S21W52 L=251)   398     30        2 (1656/ 587-1698/ 631)  14.0 W50-W54 S22-S19
     13  352  570 (S33E53 L=146)   706     80        1 ( 310/ 554- 393/ 592)   0.0 E58-E47 S33-S30
     14 1336  572 (S25W25 L=223)   247     80        0 (1303/ 556-1378/ 587)  11.3 W22-W27 S25-S23
     15 1050  478 (S33W06 L=204)   290     69        0 (1014/ 464-1080/ 498)  26.6 W03-W07 S34-S31
     16 1188  414 (S37W18 L=216)   342     85        0 (1146/ 395-1237/ 430)  23.2 W14-W21 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