Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  959 1792 (N38E32 L=279)   397     88        2 ( 913/1765-1002/1815)   0.0 E36-E27 N36-N39
      2 1160 1694 (N37E13 L=298)   488    123        1 (1104/1672-1214/1723)   0.0 E18-E08 N35-N39
      3  457 1540 (N09E51 L=260)   600    136        1 ( 427/1487- 506/1590) -24.4 E56-E44 N06-N11
      4 1633 1544 (N45W27 L=339)   384     53        0 (1601/1528-1660/1564)  26.6 W24-W30 N42-N46
      5 1653 1451 (N38W30 L=341)   558    109        0 (1622/1395-1671/1498) -56.3 W26-W31 N34-N41
      6 1594 1282 (N25W27 L=338)   412    166        1 (1529/1238-1629/1314)  -7.1 W22-W30 N22-N26
      7 1146 1235 (N08W01 L=312)   186     60        2 (1131/1198-1162/1265) -78.7 W00-W01 N05-N10
      8 1799 1121 (N22W47 L=359)   927    260        1 (1745/1063-1846/1200) -21.0 W40-W53 N18-N25
      9  869  982 (S13E07 L=304)   357    160        0 ( 819/ 924- 908/1056) -24.0 E12-E03 S15-S10
     10 1723  969 (N11W44 L=355)   635    123        2 (1676/ 946-1772/ 994)  12.5 W39-W48 N08-N13
     11  795  905 (S19E10 L=302)   308    118        0 ( 736/ 880- 844/ 926)   8.1 E13-E06 S19-S17
     12 1665  877 (N04W42 L=353)   278     82        1 (1626/ 864-1704/ 899)  38.7 W39-W44 N02-N06
     13  165  711 (S46E72 L=239)   740     79        1 ( 137/ 683- 200/ 736)   2.4 E88-E64 S47-S44
     14 1286  707 (S16W24 L=335)   272     73        1 (1250/ 691-1326/ 723)  21.8 W21-W26 S17-S14
     15  479  578 (S48E27 L=284)   332     67        0 ( 458/ 551- 506/ 607) -26.6 E29-E23 S49-S45
     16 1695  499 (S13W64 L= 15)   682     94        2 (1672/ 453-1714/ 548) -39.8 W59-W68 S15-S10
     17 1301  437 (S30W38 L=349)   537    172        2 (1235/ 384-1372/ 492)   0.0 W30-W45 S31-S28
     18 1519  380 (S26W58 L=  9)   556     86        1 (1477/ 351-1560/ 404)   0.0 W52-W64 S26-S24
     19  946  246 (S56W31 L=342)   476     68        0 ( 930/ 204- 962/ 289) -26.6 W25-W37 S58-S52
  

    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