Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  387 1549 (N16E54 L= 32)  1064    208        2 ( 340/1500- 433/1594) -15.3 E60-E49 N12-N19
      2 1743 1521 (N41W46 L=133)   509     94        0 (1719/1495-1766/1545)  45.0 W42-W49 N38-N43
      3  330 1463 (N10E56 L= 30)   389     67        1 ( 311/1428- 356/1489)  90.0 E58-E53 N07-N11
      4 1570 1416 (N28W29 L=115)  1104    219        2 (1457/1349-1649/1472) -12.1 W20-W35 N24-N29
      5 1522 1141 (N09W28 L=114)   806    171        1 (1488/1035-1548/1220) -63.4 W25-W31 N03-N14
      6 1736 1068 (N10W45 L=131)   389    102        3 (1701/1030-1770/1102) -31.0 W41-W47 N08-N11
      7  293 1056 (S15E49 L= 37)   327     62        1 ( 267/1029- 316/1084) -21.8 E51-E46 S16-S13
      8 1146  763 (S20W12 L= 99)   738    239        1 (1043/ 702-1263/ 837)  47.7 W07-W17 S24-S13
      9 1776  599 (S15W63 L=150)   509     78        1 (1753/ 574-1803/ 621)  63.4 W60-W65 S16-S12
     10 1715  509 (S21W62 L=148)   512     51        1 (1687/ 488-1740/ 533)  31.0 W58-W63 S22-S19
     11 1090  355 (S47W22 L=109)   688    125        2 (1026/ 331-1147/ 384)   5.2 W16-W27 S48-S44
     12  862  206 (S67W11 L= 98)   491     49        1 ( 834/ 181- 891/ 223)   7.1 W06-W14 S69-S64
  

    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