Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1598 1637 (N26W59 L=200)   572     51        0 (1573/1611-1628/1662)   6.3 W54-W63 N25-N27
      2  918 1341 (N22E00 L=141)   210     23        1 ( 896/1321- 939/1364) -63.4 E01-W00 N20-N23
      3  181 1268 (N32E64 L= 77)   512     57        0 ( 155/1241- 210/1288)  14.0 E68-E59 N30-N32
      4  251 1233 (N28E53 L= 88)   524     57        2 ( 225/1206- 270/1267)  63.4 E55-E50 N26-N30
      5 1612  929 (S15W35 L=176)   389    103        0 (1577/ 886-1638/ 970)  38.7 W31-W37 S16-S12
      6 1796  900 (S21W51 L=192)   988    198        1 (1713/ 842-1863/ 962)  13.2 W42-W59 S22-S17
      7  484  618 (S12E43 L= 98)  1072    160        1 ( 456/ 526- 517/ 698)  48.0 E47-E38 S17-S07
      8 1705  649 (S35W43 L=184)   332     69        1 (1679/ 627-1728/ 671) -39.8 W39-W46 S37-S32
      9 1451  567 (S35W20 L=161)   789    226        2 (1360/ 535-1565/ 612) -29.4 W13-W29 S39-S29
     10 1597  554 (S39W33 L=174)   369     87        1 (1571/ 528-1623/ 600)  45.0 W30-W35 S40-S36
     11 1222  471 (S36W01 L=142)   846    216        1 (1138/ 407-1298/ 538)  22.6 E05-W07 S39-S32
     12 1069  482 (S32E09 L=132)   652    124        0 (1042/ 429-1101/ 533)  74.1 E11-E05 S35-S28
  

    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