Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  758 1652 (N54W00 L=356)   505     72        1 ( 719/1624- 788/1679) -38.7 E03-W03 N50-N55
      2 1576 1609 (N24W59 L= 54)  2095    183        5 (1510/1559-1657/1660) -71.6 W51-W67 N18-N27
      3 1634 1526 (N18W58 L= 53)   538     53        2 (1609/1500-1665/1551)   0.0 W54-W61 N15-N19
      4 1425 1468 (N21W38 L= 34)   945    253        1 (1303/1410-1527/1511) -45.0 W31-W42 N15-N26
      5 1290 1307 (N16W23 L= 18)   761    262        5 (1243/1213-1352/1365)  37.9 W18-W27 N11-N19
      6 1578 1141 (S01W36 L= 31)  1059    200        3 (1517/1064-1618/1222)  33.7 W30-W41 S04-N02
      7  736 1154 (N21E14 L=341)   231     61        0 ( 716/1122- 756/1188)  45.0 E15-E12 N19-N22
      8  196 1081 (N27E62 L=293)   508     75        2 ( 164/1048- 228/1108) -24.0 E67-E57 N25-N29
      9 1707  901 (S19W40 L= 36)   292     74        0 (1664/ 886-1735/ 913) -33.7 W36-W42 S20-S16
     10 1537  689 (S26W24 L= 19)   792    192        1 (1477/ 646-1608/ 757) -51.3 W20-W29 S30-S20
     11 1642  715 (S28W33 L= 28)   339     79        1 (1614/ 687-1669/ 754)  36.9 W29-W35 S29-S26
     12 1814  691 (S37W52 L= 48)   869     85        1 (1768/ 663-1855/ 723)  -4.4 W45-W59 S38-S34
     13 1728  592 (S41W42 L= 37)   509     51        0 (1693/ 576-1761/ 606) -18.4 W37-W46 S42-S38
     14 1660  536 (S42W34 L= 30)   363     32        0 (1640/ 515-1682/ 556)  26.6 W32-W36 S43-S40
     15 1571  492 (S42W25 L= 20)   417     75        0 (1542/ 473-1608/ 518) -45.0 W22-W28 S44-S38
     16  977  342 (S34E24 L=332)   601     92        0 ( 940/ 317-1024/ 378) -59.0 E25-E21 S36-S30
     17  546  326 (S25E57 L=298)   862    122        1 ( 522/ 287- 585/ 373)  22.6 E63-E50 S28-S22
     18  621  298 (S29E53 L=303)   833     60        1 ( 597/ 264- 649/ 335)  63.4 E56-E49 S31-S25
  

    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