Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1135 1797 (N51W43 L=286)   330     56        1 (1107/1777-1156/1813)  76.0 W40-W46 N49-N53
      2  697 1440 (N40E11 L=232)   689    170        2 ( 626/1416- 773/1476) -18.4 E17-E05 N36-N42
      3 1268 1441 (N24W27 L=270)   909    315        0 (1146/1417-1410/1474) -18.4 W18-W36 N20-N27
      4  576 1425 (N42E22 L=222)   243     68        0 ( 551/1396- 596/1451) -53.1 E23-E20 N40-N44
      5 1876 1211 (S07W67 L=310)   531     61        1 (1849/1181-1895/1231) -53.1 W63-W70 S08-S04
      6  764 1173 (N22E12 L=231)   604    232        0 ( 674/1112- 854/1223) -45.0 E17-E07 N16-N26
      7  342 1013 (N21E46 L=197)   633    143        2 ( 289/ 974- 397/1064) -41.2 E51-E41 N17-N24
      8 1909  931 (S24W66 L=309)  1415    180        2 (1868/ 865-1939/1022)  28.3 W58-W73 S28-S19
      9  686  971 (N12E21 L=222)   217    126        0 ( 628/ 960- 751/ 979) -29.7 E24-E17 N09-N13
     10 1446  696 (S23W18 L=261)   450    190        1 (1371/ 638-1504/ 762)  20.0 W11-W22 S24-S20
     11  571  686 (S03E36 L=208)   583     89        1 ( 528/ 648- 606/ 723)  26.6 E39-E33 S04-N00
     12 1750  604 (S41W45 L=289)   541     53        1 (1717/ 576-1779/ 620) -21.8 W40-W50 S43-S37
     13  564  502 (S14E43 L=200)   671    143        1 ( 522/ 444- 620/ 558) -77.5 E44-E41 S18-S09
     14 1678  501 (S46W38 L=281)  1480    150        1 (1623/ 458-1741/ 551)  13.0 W30-W46 S49-S42
     15 1494  375 (S49W17 L=261)   612    106        2 (1461/ 343-1537/ 418)  29.7 W13-W22 S51-S46
     16 1407  332 (S50W08 L=252)   347     40        0 (1394/ 310-1425/ 359) -90.0 W06-W09 S51-S46
     17 1111  295 (S42E17 L=227)   687    103        2 (1074/ 267-1148/ 329)  -9.5 E19-E13 S44-S38
     18  627  282 (S30E54 L=190)   524     63        1 ( 590/ 267- 660/ 299) -21.8 E58-E50 S31-S29
  

    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